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ABSTRACT 

A  theoretical  study  of  the  gas  dynamic  interaction  between  a 

hypersonic  plume  and  the  opposed  hypersonic  external  stream  is 

presented.  Steady,  axisymmetric,  inviscid,  perfect  gas  flow  is 

postulated  for  both  the  bow  and  far  field  regions.  Limiting 

forms  of  solutions  are  obtained  for  the  bow  region  by  application 
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ambient  air  flow.  Through  asymptotic  expansions  and  their  matching, 
it  is  found  that  six  regions  are  required  to  adequately  describe 
the  bow  region.  For  the  far  field  region,  the  hypersonic  small- 
disturbance  form  of  the  Nevton-Busemann  approximation  (i.e.,  e  ->0, 
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expansions  and  their  matching,  it  is  found  that  "entropy  wake" 
solutions  are  required  to  adequately  describe  the  exhaust  flow  and 
the  air  flow  near  the  contact  surface.  Analytical  solutions  are 
obtained  which  (i)  define  scaling  parameters  for  the  bow  and  far 
field  flow;  (ii)  estimate  the  accuracy  of  the  Newtonian  impact 
theory  in  predicting  bow  region  geometry  and  properties;  (iii)  estab¬ 
lish  the  variation  of  bow  and  far  field  properties  with  variation 
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INTRODUCTION 


A  theoretical  study  of  the  interaction  of  a  highly  underex¬ 
panded  jet  issuing  into  a  high  Mach  number  free  stream  that 

* 

is  opposed  to  the  jet’s  expansion  along  its  aocis  is  present¬ 
ed.  At  some  distance,  characterized  by  the  length  R.aio 
the  jet  and  free  stream  gases  interact.  The  jet  flow  is  there 
by  deflected  downstream  by  the  external  flow.  The  particular 
case  of  interest  to  be  examined  in  this  study  concerns  flows 
where  the  typical  dimension  of  the  body  D  is  much  less 
than  RsCo  (see  schematic  (1) ) . 


Schematic  (1) 


For  this  case,  analytical  models  of  the  upstream  and  downstream 
gas  dynamic  interaction  processes  are  formulated  herein.  The 
following  approximations  are  postulated: ’ 

1.  Steady  Flow 

2.  Axisymmetric  Flow 
3«  Continuum 

\  4.  Inviscid 

5.  Thermally  and  Calorically  Perfect  Gas 

In  further  studies,  assumption  4  can  be  partially  relaxed  by 
applying  boundary  layer  concepts  along  the  shear  surfaces  for 
high  Reynolds'  number  flows.  Also,  assumption  5  can  be  re¬ 
laxed  quite  easily  by  applying  Mollier  charts  or  similar  models 
for  equilibrium  gas  dynamic  properties. 


A  general  schematic  of  the  flow  structure  under  these  assump¬ 
tions  is  given  in  figure  1,  where  three  distinct  flow  regions 
are  discerned.  The  major  characteristics  of  these  regions 
are  as  follows: 

Bow  Region  ( 

This  region  is  characterized  by  nearly  normal  interior  and 
exterior  shocks  and  relatively  thin  layers.  The  flow  within 
the  layers  is  then  necessarily  subsonic  and  of  high  density. 
Corner  Region 

The  layers  thicken  markedly  and  turn  in  a  downstream  direc- 

t 

tion.  Both  interior  and  exterior  layers  undergo  transonic 
expansions  about  the  relatively  motionless  core  region  sur¬ 
rounding  the  plume.  The  development  of  the  plume  ahead  of 
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the  internal  shock  resembles  very  closely  that  of  a  plume  ex¬ 
panding  into  a  quiescent  ambient. 


Far  Field  Region 


The  exterior  shocrc  and  interface  are  at  relatively  small  angles 
with  the  freestream.  The  flow  in  all  layers  is  supersonic  and 
nearly  parallel  to  the  axis. 


Flow  systems  of  these  types  are  of  considerable  interest  for  the 
practical  applications  of  force  vector  control  and  attendant 
surface  thermal  protection  of  re-entry  vehicles.  For  force  vec¬ 
tor  co'.  rol,  the  altered  pressure  distribution  on  the  body,  as 
well  as  the  jet  thrust,  has  to  be  considered  in  calculating  the 
total  force  on  the  body.  For  thermal  protection,  although  the 
hot  external  gases  are  blown  free  from  the  nose  by  the  jet,  it 
is  still  necessary  to  consider  the  heat  flux  due  to  the  poss¬ 
ible  reattachment  of  this  separated  gas  flow  on  the  body. 

In  addition  to  the  effects  already  mentioned,  which  are  mani¬ 
fested  near  the  body,  we  must  also  be  concerned  with  the  flow 
field  far  downstream  of  the  body.  In  this  wake-like  region, 
the  mixing  of  the  jet  gases  may  affect  the  chemical  processes 
occurring  in  the  wake  to  such  an  extent  that  observable  prop¬ 
erties  and  hence  detection  or  communication  may  be  greatly  fa¬ 
cilitated  or  decreased.  It  is  apparent  that  if  such  a  sys¬ 
tem  were  to  be  properly  utilized,  a  detailed  understanding  of 
the  physical  and  chemical  processes  occurring  within  the  dis¬ 
turbed  flow  field  must  be  obtained. 


Previous  theoretical  and  experimental  studies  have  added  con¬ 
siderably  to  our  basic  knowledge  in  this  area.  The  first  stud- 
12  3 

ies  *  *  were  essentially  experimental  in  nature,  applying  only 
simplified  modeling  to  correlate  the  data.  A  significant  result 
of  these  early  investigations  was  that  two  modes  of  inter¬ 
action  were  observed  to  exist.  The  modes  were  dependent  upon 
the  jet  exit  Mach  number  Mj  ,  free  stream  Mach  number  Moot  , 
ratio  of  jet  stagnation  pressure  to  free  stream  dynamic  pres¬ 
sure  TT  ,  body  size  D  ,  and  shape.  One  mode  is  characterized 
by  large  interaction  distances  with  unsteady  shocks  and  bound¬ 
ing  surfaces.  This  flow  regime  is  designated  as  the  unstable 
case . 

The  second  mode  of  interaction  is  characterized  by  a  relatively 
short  interaction  length  and  steady  strong  shocks.  This  flow 
regime  is  the  case  investigated  herein.  The  mechanism  for 
transition  from  the  steady  to  unsteady  flow  was  postulated  by 
Finley  to  be  a  result  of  the  development,  for  low  Tf  and 

or  high  Mj  ,  of  a  multiple  cell  structure  for  the  un¬ 
disturbed  plume  before  the'  interaction  region  is  reached.  In 
view  of  the  subsonic  flow  existing  behind  the  Mach  disc  sepa¬ 
rating  the  cells,  we  have  the  possibility  of  upstream  influ¬ 
ence  from  the  surrounding  gases  affecting  the  interaction  re¬ 
gion  in  a  possibly  unsteady  manner. 
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The  boundary  between  stable  and  unstable  flow,1’2’ 3for  experi 
ments  where  the  effect  of  body  size  on  the  flow  is  negligible 
is  given  by: 


Mooe 

Mj 

Tf 

7.1 

1.0 

2 

7.1 

4.0 

8 

7.1 

4. 35 

15 

7.1 

5.30 

20 

2.71 

3.10 

4 

2.71 

3.90 

8 

Thus,  to  achieve  stable  flow,  the  jet  stagnation  pressure  must 
be  much  larger  than  the  free  stream  dynamic  pressure  (i.e.  ,Tf 
>V  1  )  for  moderate  values  of  Mj  .  These  are  the  same  con¬ 
ditions  under  which  R*,.',,  becomes  much  larger  than  D  (for  D 
approximately  equal  to  jet  exit  diameter);  therefore,  the  case 
studied  herein  will  always  be  of  the  stable  type. 

It  should  be  noted  that  one  investigator  2  observed  a  region 
for  large  values  of  TT  where  the  flow  became  unstable  and 
continued  to  be  so  for  all  higher  values  of  TT  .  This  result 
has  not  been  duplicated  by  other  investigators,  and  indica¬ 
tions  of  stable  interactions  ^  have  been  obtained  from  other 
experiments  at  conditions  that  fall  within  the  region  of  un¬ 
stable  flow  found  in  Reference  2. 

In  addition,  it  has  been  observed  that  the  interaction 

length  for  the  stable  condition  is  a  function  of  TY,  Vj 

and  that  the  contact  surface  separating  the 
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external  stream  from  the  jet  gases  is  almost  spherical  in  nat¬ 
ure.  These  observations  have  been  utilized  in  constructing 
simplified  models  to  describe  the  flow  field  associated  with 
stable  interactions. 

Theoretical  analyses  of  this  jet  interaction  problem  have  ap¬ 
peared  only  recently.  In  1969,  Laurmann  6, utilizing  the 
Newtonian  impact  theory,  calculated  quantitative  results  for 
the  upstream  region.  This  paper  presented  general  features  of 
the  upstream  bow  interaction,  but  did  not  correctly  treat  such 
areas  as  the  corner  region  and  the  far  field  development  of 
the  flow.  In  1971,  a  time-dependent  numerical  technique  ^  was 
employed  by  Rudman  and  Vaglio-Laurin  to  calculate  detailed 
quantitative  dynamic  properties  of  the  upstream  flow  interac¬ 
tion.  In  the  calculations,  the  jet  plume  boundary  layer  in¬ 
tersection  with  the  shock  layers  was  assumed  to  be  a  point  in¬ 
teraction  and  the  core  pressure  was  assumed  to  be  that  of  the 
freestream.  The  downstream  flow  was  not  considered  in  the 
above  report  and  only  the  bo w  region  was  evaluated. 

The  present  study  was  undertaken  for  the  primary  task  of  cor¬ 
recting  the  errors  inherent  in  a  Newtonian  impact  analysis,  to 
estimate  the  effect  of  the  corner  region  on  the  total  flow 
field,  and  to  calculate  the  far  field  interaction. 

An  analytic  approach  was  considered  since  such  a  method  pro¬ 
vides  for  the'  greatest  insight  into  the  physical  processes  and 
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Also  defines  the  accuracy  associated  with  various  approxima¬ 
tions.  Finally,  this  study  also  acts  as  a  model  for  extend¬ 
ing  calculations -to  more  complicated  configurations;  e.g., 
where  the  body  influences  the  flow  or  where  the  jet  is  at  an 
angle  of  incidence  to  the  free  stream. 

Vie  must  now  determine  the  operating  conditions  under  which  our 
simplifying  assumptions  are  satisfied.  For  a  missile  or  ve¬ 
hicle  re-entering  the  earth's  atmosphere,  we  have  the  follow¬ 
ing  range  of  conditions  from  approximately  400,000  ft.  alti¬ 
tude  down  to  approximately  50,000  ft.  altitude. 

Range  of  Typical  Operating  Conditions  for  a  Re-Entry  Vehicle 

Uoo*  (ft/sec)  15,000  23,000 

Mco«  10  25 

poo*  (lbf/ft2)  6.92xlO“5  23.27 

Over  this  range  of  operating  conditions,  we  generally  find 
that  rarefaction,  viscous,  chemical  or  physical  effects  may  be 
significant.  This  would  violate  our  assumptions  and  therefore 
we  must  examine  the  relevant  parameters  and  their  numerical 
values  in  order  to  ascertain  when  our  assumptions  are  fully 
satisfied. 

The  first  limitation  on  the  parameters  involves  the  assump¬ 
tion  of  steady  flow  and  requires  that  1 ,  M*,**  1  j  Mj  c  0(l) 

The  reasons  for  this  have  been  discussed  previously.  A  second 
limitation  relates  to  rarefaction  effects  occurring  within  the 
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flowfield.  For  the  upstream  region,  we  have  three  areas  where 
rarefaction  effects  may  invalidate  the  application  of  the  ana¬ 
lysis.  The  first  area  is  the  plume  core  itself  where  the  large 
degree  of  expansion  may  take  the  free  plume  flow  into  the  tran¬ 
sitional  regime.  From  an  analysis  of  a  steady  spherical  source 

d 

flow  expanding  into  a  vacuum  Kamel  and  Y/illis  derived  an  ex¬ 
pression  for  the  distance  beyond  which  the  source  gas  becomes 
collisionless.  This  is  given  by 

R4  (Re*)^ 

where  is  the  jet  nozzle  throat  radius  and  Re*  is  the  Rey¬ 
nolds'  number  evaluated  at  the  nozzle  throat.  This  value  is 
only  a  function  of  the  nozzle  throat  conditions  and  hence  is 
a  constant  for  a  particular  set  of  chamber  conditions.  This 
can  also  be  simply  related  to  the  thrust  of  our  jet 


where  is  the  specific  impulse,  T  is  the  thrust  and  R>,'  is 
the  jet  chamber  pressure.  In  order  that  the  plume  flow  for 
our  problem  be  a  continuum,  the  interaction  distance,  Rsta  , 
for  the  interior  shock  must  be  much  less  than  the  distance  R-t- 
The  interaction  distance  Rst'0  is  determined  by  the  condi¬ 
tion  of  equal  stagnation  pressures  along  the  axis.  Therefore, 


This  expression  is  a  result  of  assuming  a  source  like  behavior 
for  the  undisturbed  plume  properties.  The  occurence  of  the 
various  quantities  in  the  above  expression  will  be  detailed 
in  a  later  section  relating  to  undisturbed  plume  properties. 
The  condition  for  continuum  flow  for  the  plume  expansion  is 
then 

B±-  =.  (Re*)  IjL  ^  1 

Rsio  ^St’6 

In  terms  of  the  thrust  of  the  jet,  the  above  becomes 


Following  Reference  (9)  we  take  the  numerical  value  of  —  10 

Rst’o 

to  be  safely  within  the  continuum  region.  This  bound  is  in¬ 
dicated  in  Figure  2  for  the  representative  conditions  of 

Meoe  -  10  J  =  50  atm.  >yu*  =  zux  io“6  Ifcf  s.c/ft*,  U*-3000fpS 
tfi.  =  1«2.5  )  Xsp  —  300  **o 

We  now  direct  our  attention  to  the  internal  and  external  shock 
layers  where  rarefaction  effects  can  result  in  the  thickening 
of  the  shocks  and  shear  layers  to  such  an  extent  that  they 
strongly  influence  the  development  of  the  layers.  In  a  paper 
by  Bush10  it  was  illustrated  that  the  thickening  of  the  shocks 
and  growth  of  the  viscous  shear  layers  are  related  for  the 
shock  layers;  and,  therefore,  we  need  only  determine  the  thick¬ 
ness  of  the  shock  with  respect  to  the  layer  thickness  to  de- 


jy.  pi  mi': 


termine  where  inviscid  continuum  theory  is  valid.  It  has  been 
11 

shown  that  a  continuum  inviscid  description  of  the  flow  is 
appropriate  when' 

>  s^/stl  -  15 

where  in  the  above  £  denotes  thickness,  and  X  denotes  layer, 

S  snock  and  ®  ,  l  ,  external  and  internal  layers  respect¬ 
ively.  From  Adams  and  Probstein  12  we  have  for  the  thickness 
of  the  shock  for  either  external  or  internal  flow 


Ss  -  V>*/  C# 

Where  C*  and  V>*  respectively  are  the  speed  of  sound  find 
cosity  evaluated  at  M  =  1  . 


vis- 


After  introducing  the  interaction  length  Rs.’c  as  a  signifi¬ 
cant  length,  this  then  reduces  to  ^ 


4s-  -  ~  R 

Rsio 


oo 


Where  UJ  is  the  exponent  for  temperature  in  the  power  law 
viscosity  model,  R®  is  the  Reynolds’  number  arid  Mco  is  the 
Mach  number,  based  on  free  stream  properties.  We  now  intro¬ 
duce  the  first  order  thicknesses  of  the  shock  layers,  which 
will  be  derived  in  detail  in  a  later  section. 


Rs 


_  €, 


i© 


So  s-  Vi 

Rsi'o 
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Uhere  6«,£io  are  the  density  ratios  across  the  normal  part  of 
the  shocks  along  the  axis  of  symmetry. 

Therefore , 

4^  —  £“  R®«  CZ._U*taL  Roo*' 

Ss«  ^Uoc.  dsc  ZZZtt 

where  Ra>  is  based  on  the  interaction  distance  Rsi0  .  Using 
the  source  model  for  the  plume  properties  and  the  equality  of 
stagnation  pressures  to  determine  Rs.c  ,  we  can  reduce  the 
above  ratios  to  functions  of  jet  thrust,  ~T  ,  jet  throat  prop¬ 
erties  and  external  flow  gas  properties.  The  results  of  this 
are  • 

jL&di  =  -£c  Floo*  C~T  pooe)  _ _  o 

u“e  r«pVx  <£•  VH  (  j£~c  05  "ggjM*  fd<p)  ^  ” 

hk  =  €  *  *  *  M*fc  *<f  -  , 

dSc  I5pVi 

Using  values  stated  previously  which  are  representative  of 
typical  flight  conditions,  the  above  can  be  plotted  as  a  func¬ 
tion  of  altitude  for  a  given  thrust,  above  which  rarefaction - 
effects  have  to  be  accounted  for.  This  has  been  carried  out 
in  Figure  2  and  it  can  be  seen  that  for  jets  of  10  to  10^ 
pound  force  thrust,  there  is  a  large  region  of  flow  where  the 
continuum  inviscid  analysis  will  apply. 

V/e  must  no w  examine  the  region  where  perfect  gas  behavior  or 
equilibrium  gas  behavior  will  apply.  For  the  external  shock 
layer,  there  are  many  analyses  and  data  available  which  indi¬ 
cate  when  nonequilibrium  effects  in  such  general  flow  fields 


become  important.  For  the  plume  and  internal  flow,  this  is  not 
the  case.  The  reason  for  this  is  that  rocket  exhaust  chemi¬ 
cal  composition, 'especially  those  of  the  solid  propellant  va¬ 
riety,  varies  greatly  from  application  to  application  and  in 
each  case  the  bounds  of  equilibrium  are  different.  Therefore, 
we  cannot,  define  a  single  general  condition  when  nonequili¬ 
brium  flow  occurs.  Therefore,  we  will  set  down  the  bound  for 
the  external  flow  as  the  limit  for  the  combined  layers  above 
which  the  flow  cannot  be  considered  to  be  in  equilibrium. 

The  external  layer  nonequilibrium  bound  is  adapted  from  a 
paper  by  Cheng14  and  is  based  on  the  relaxation  processes  be¬ 
hind  a  normal  shock.  When  the  distance  for  relaxation  to  equi¬ 
librium  conditions  behind  the  shock  becomes  of  the  same  order 
of  magnitude  as  the  layer  thickness,  then  nonequilibrium  ef¬ 
fects  must  be  taken  into  account.  Cheng's  results  gave  a  sin¬ 
gle  point  at  high  altitude,  which  when  combined  with  binary 
scaling  led  to  the  result  in  Figure  2,  above  which  nonequili-  • 
brium  effects  must  be  considered.  Binary  scaling  derives  from 
the  fact  that  at  high  altitude  the  probability  for  chemical 
reactions  to  occur  by  three-body  collisions  is  much  less  than 
that  for  two-body  reactions.  Under  these  conditions,  it  can 
be  shown  that  if  density  and  field  dimension  are  held  constant, 
then  the  degree  of  nonequilibrium  will  be  the  same  in  each 
case . 
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Since  the  body  dimension  ^s4‘4is  related  to  jet  thrust  and  den¬ 
sity  pae  is  related  to  altitude,  we  can  then  establish  a 
relationship  between  thrust  and  altitude,  which  will  maintain 
the  product  of  density  and  dimension  constant,  giving  the 
same  degree  of  nonequilibrium.  The  binary  scaling  principle 
only  applies  to  high  altitude,  for  moderate  to  low  altitude 
collisions  should  be  prevalent  enough  that  equilibrium  flow 
is  maintained. 

Referring  to  Figure  2,  it  is  shown  that  for  typical  re-entry 
conditions  the  extent  of  the  equilibrium,  inviscid,  continuum 
flow  is  about  one  half  of  the  complete  continuum,  inviscid 
region;  however,  there  is  still  a  considerable  region  where 
all  the  assumptions  stated  previously  are  fully  satisfied. 

Now  that  we  have  determined  our  bounds  of  validity  for  an  ana¬ 
lysis  based  on  several  assumptions,  we  can  proceed  to  outline 
the  analysis  in  detail. 
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JET  EXHAUST  MODEL 

The  undisturbed  plume  which  interacts  with  the  external  flow 
to  form  the  bow  layer  develops  in  an  environment  (the  core 
region)  where  there  is  little  fluid  motion  and  hence  nearly 
constant  pressure.  Under  these  circumstances,  the  undisturbed 
plume  is  analogous  to  the  case  of  a  jet  exhausting  into  a 
quiescent  ambient  which  is  at  the  core  pressure.  V/e  can  then 
employ  analyses  developed  for  plumes  exhausting  into  quies¬ 
cent  ambients  to  our  retro  plume  case.  To  be  consistent  with 
the  analyses  for  the  downstream  and  bow  region,  vie  will  uti¬ 
lize  an  analytical  approach  for  the  undisturbed  plume.  There 
are  many  models  available  for  describing  the  isentropic  ex¬ 
pansion  of  a  jet  far  away  from  the  nozzle  exit.  Many  of  them 
have  the  form 


(1) 


Other  properties  can  be  determined  from  this  expression  by  use 
of  the  isentropic  flow  relations  (No’s.  30.2-30.3  from  the  Bow 
Analysis) . 


*In  this  reference,  an  expression  for  Mach  number*  is  given 
vihich  when  substituted  into  the  isentropic  relationships  for 
high  Mach  number  gives  a  result  identical  to  (1). 
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In  (1)  <p  is  the  angle  measured  away  from  the  axis  of  sym¬ 
metry,  Rv  is  radius  measured  from  the  nozzle  exit  and  n 
and  C  are  constants  for  a  given  nozzle  and  jet  gas  compo¬ 


sition.  Various  forms  have  been  proposed  for  H  and  C  as 
functions  of  plume  specific  heat  ratio  ,  exit  Mach  num¬ 
ber  Mj  ,  and  nozzle  exit  angle  ©n  .  Using  the  form  C  = 
7T/20OO  , where  0oo  is  the  sum  of  the  Prandtl-Meyer  limiting 
turning  angle  into  a  vacuum  and  ©n5gives  the  correct  theo¬ 
retical  result  of  zero  density  or  infinite  Mach  number  when 

9=  ©  oo  .  For  this  value  of  C  ,  we  can  extract  from  a 

15 

Prandtl-Meyer  analysis  near  the  nozzle  lip  the  value  f or n. 


n  » 


*i-  l 


(2) 


This  result  can  also  be  obtained  from  application  of  the  small 

1C 

disturbance  theory  to  hypersonic  jets  expanding  into  a  vacuum  . 
As  will  be  shown  in  the  bow  layer  analysis  section,  these  va¬ 
lues  for  H  and  C  give  the  physically  realistic  resul't  in  the 
limit  as  Gto  O  that  the  density  and  other  flow  vari¬ 

ables  remain  of  the  same  order  as  we  go  away  from  the  axis. 

For  any  other  choice  of  C  and  H  as  a  function  of  <$j  the 
density  would  be  finite  on  the  axis  and  either  zero  or  infi¬ 
nite  away  from  the  axis.  Even  though  the  form  of  n  as  a 
function  of  Kj  is  correct,  its  numerical  value  for  a  given 
is  still  not  clear.  This  is  because  comparisons  with  numerical 
results  give  different  values  for  H  ,  varying  between  the  value 


» 


m 


given  in  (2)  and  one  half  that  value  The  difficulty  is 

that  the  form  of  equation  (1)  is  not  general  enough  to  be 
valid  over  a  full  range  of  R  and  (f  .  This  point  is  illus¬ 
trated  in  an  article  by  Boynton  15  where  he  shows  for 
the  proportionality  constant  is  approximately  one  and  for 
9>  ,lf  the  numerical  results  lie  closer  to  the  curve  for  a 
proportionality  factor  of  two.  Since  the  bow  shock  layer  and 
downstream  layer  depend  critically  on  the  undisturbed  plume 
flow,  (as  will  be  shown  in  the  bow  and  downstream  analyses 
sections)  then  it  was  felt  that  a  review  of  numerically  cal- 
*  culated  undisturbed  plumes  to  give  "best"  values  for  H  was 


justified.  The  expected  value  would  probably  lie  somewhere 

between  the  theoretical  value  . . -  and  - . . — 

-  3L  cUlQ  -  l 

Three  sources  of  numerical  plume  calculation  data  were  chosen 
for  examination  to  determine  "best"  values  for  n  .  In  re¬ 


ference  (20)  numerical  calculation  of  three  plumes  from  super¬ 
sonic  nozzles  expanding  into  a  quiescent  ambient  were  presented. 
Two  of  the  nozzles  (  -  1.15,1.18)  are  representative  of 

those  used  for  launch  vehicles,  and  the  third  (  Xj  *  1.2*0  ±s 
representative  of  spacecraft  nozzles.  The  method  of  charac¬ 
teristics  was  used  in  this  study  to  determine,  for  constant 
J  ,  Mach  number  contours,  and  density  and  temperature  dis¬ 
tributions  along  and  perpendicular  to  the  axis  of  symmetry. 

The  second  reference  16  contained  the  calculation  of  the  plume 
resulting  from  the  expansion  of  air  through  a  sonic  orifice. 


Mach  number  distributions  perpendicular  to  the  axis  of  sym- 


metry  were  given.  The  third  source  used  15  contains  results 
of  a  numerical  calculation,  using  a  finite  difference  tech¬ 
nique,  of  the  flow  from  a  nozzle  representative  of  a  launch 
vehicle.  Density  is  given  as  a  function  of  angle  from  the 
axis  in  this  work.  To  determine n from  the  above  sources  for 
points  away  from  the  axis,  we  simply  use  the  isentropic  re¬ 
lationships  and  the  assumed  distribution  (1)  to  determine 
for  a  known  value  of  the  constant  ©«,  and  jo  at  a  given  <p 
and  R  .  However,  near  the  axis  this  procedure  does  not 
yield  good  results  since  large  changes  in  n  produce  little 
measured  change  in  ^  or  h  in  this  region.  In  this  case, 
we  then  apply  a  derivative  of  the  density  distribution  to  ob¬ 
tain  a  better  relationship  for  determining  n  .  For  the  first 
and  second  references  where  properties  are  given  as  functions 

of  T  ,  the  distance  perpendicular  to  the  axis,  the  expres¬ 
sion 


establishes  n  from  the  density  gradient  perpendicular  to  the 
axis.  For  the  third  reference,  we  use 


(4) 
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Values  of  n  »  the  ratio  of  n  calculated  from  3  and  4,  to 
2/(*j'-.L)  are  plotted  tin  figure  (3)  as  a 'function  of  $>/©«, and 
X/rj  .  From  this  graph,  it  is  evident  that  n  is  not  a  con¬ 
stant  but  varies  both  as  a  function  of  y/Qec  and  X/r :  .  If 
use  is  to  be  made  of  (1)  in  calculating  retro  plumes,  then 
for  each  change  in  X/ r j  a  new  value  of  ft  must  be  determined 
from  figure  (3). 

The  variation  of  ft  with  ^>/0ttis  not  as  critical  as  the  axial 
variation  since  the  bulk  of  the  mass  flow,  which  determines 
the  plume  external  stream  interaction,  lies  near  the  axis  and 
therefore  the  behavior  for  5>/©00>  *5  is  not  significant  to  the 
interaction.  Also,  for  ty/Qta  near  zero  the  density  is  near 
one  for  any  value  of  ft  .  Noting  this,  we  then  take  the  value 
of  n  found  at  g)/©*,-. 2.5  to  be  the  "best"  value  for  a  given 
X/Tj  .  Characteristic  retro  plume  calculations  have  been 
carried  out  with  values  of  ft  determined  from  figure  (3).  Re¬ 
sults  of  these  calculations  will  be  discussed  in  the  sections 
on  bow  and  downstream  analyses. 


ANALYSIS  OF  BOW  REGION 


For  the  problem  studied  here  Ma>.  ,  Moo 1  and  ZC® 

The  flow  within  the  internal  and  external  shock  layers  are 

determined  by  a  perturbation  theory  for  strong  shock  waves 

1 

and  small  layer  thicknesses.  The  procedure  involves  a  limit* 
of  the  governing  equations  and  boundary  conditions  when 
€e,€Le,-+  O  Hooej  Mooi0 ->  00  and  Moof€«  ,  0(j.)where 

€e)  are  the  density  ratios  across  the  shocks  at  the 

axis.  To  elucidate  some  detail  of  the  flows  within  the  la¬ 


yers,  the  coordinates  normal  to  the  shocks  are  expanded  in 
powers  of  €e  ,  € La  where  Jt ,  m  are  to  be  determined. 

The  orders  of  magnitude  of  the  various  flow  quantities  in  terms 
of  ^6  and  Mooc^-e  j.Moo^€:L©  are  derived  from  the 


requirement  that  the  shock  relations  and  flow  equations  yield 
a  nontrivial  system  which  includes  all  the  physical  effects 
of  interest  as  6e,€t'0  ->  O  .  The  coordinate  systems  are 


shock  oriented  and  are  illustrated  in  schematic  (2) 


Schematic  (2) 

This  limit  will  be  denoted  the  N— B  limit  after  Newton— Busemann 


Upon  application  of  the,  N-B  limit  mentioned  above it  is  found 
that  the  flow  field  divides  into  several  regions  within  which 
unique  asymptotic  expansions  for  the  dependent  variables  must 
be  determined.  It  will  be  shown  in  the  following  analysis 
that  four  regions  characterized  by  two  major  effects  are  re  - 
quired  to  describe  adequately  the  external  flow  and  two  re  - 
gions  characterized  by  two  major  effects  are  needed  for  the  in¬ 
ternal  flow.  In  both  the  internal  and  external  layers,  a  re¬ 
gion  characterized  by  constant  density  to  first  order  is  formed 
near  the  axis  of  symmetry.  Away  from  the  axis,  the  density 
and  other  properties  vary  both  along  and  across  the  layers. 

For  the  external  layer,  a  region  characterized  by  non-constant 
stream  velocity  is  formed  near  the  contact  surface  in  both  the 
near  and  away  from  axis  regions.  This  contrasts  with  the  char¬ 
acter  of  the  region  near  the  shock  where  the  velocity  is  con¬ 
stant  to  first  order. 

The  regions  are  numbered  one  through  six  in  schematic  (2)*. 

Once  the  flow  properties  are  established  in  each  ofi  these  re¬ 
gions,  then  matching  between  the  expanded  properties  will  be 
demonstrated  to  show  their  consistency  and  where  necessary, 
composite  expansions  will  be  formulated.  The  final  step  in  the 
solution  to  the  bow  region  problem  will  consist  of  the  numeri¬ 
cal  matching  of  contact  surface  pressure  and  position  (i.e., 
flow  deflection)  between  the  external  and  internal  shock  layer 
flows. 

*  These  regions,  of  course,  are  not  drawn  to  scale. 


To  achieve  this  result:  (i)  initial  radii  of  curvatures  and 
positions  of  internal  and  external  shocks  will  be  found  from 
a  scheme  which  uses  solutions  obtained  for  regions  1,3>  and 
5  (this  scheme  will  be  outlined  in  appendix  (A));  (ii)  shocks 
will  be  extended  away  from  the  axis  utilizing  these  radii  of 
curvature; (iii)  external  and  internal  contact  surface  pres¬ 
sure  and  position  will  then  be  calculated,  using  solutions 
from  regions  2,  A  and  6.  The  pressures  and  positions  are  then 
compared  to  determine  if  numerical  matching  is  achieved;  (iv) 
if  equality  to  a  certain  tolerance  is  not  found,  then  the  ra¬ 
dii  of  curvatures  are  iterated  until  matching  is  achieved. 
Having  accomplished  matching,  the  procedure  starting  from 
step  (ii)  will  be  repeated  using  the  newly  found  radii  of  cur¬ 
vatures  as  the  initial  values. 

A  brief  outline  of  the  regions  considered  and  their  major 
characteristics  is  now  given. 

For  region  two,  which  is  of  0(G«)  in  thickness  and  0(1)  in 

length,  the  N— B  limit  gives  the  familiar  hypersonic  blunt  body 

result,  which  has  been  investigated  in  whole  or  in  part  by 
10  21  22 

many  authors  f  '  .  The  orders  of  magnitude  of  the  flow 

quantities  in  this  region  are  the  same  as  those  directly  be¬ 
hind  the  shock.  A  major  characteristic  of  the  solution  found 
is  that  the  velocity  along  streamlines  is  constant  in  the 
first  order  approximation.  This  result,  adequate  for  points 
near  the  shock,  has  been  shown  to  be  incorrect 


near  the 


s 
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contact  surface  in  region  4;  and,  therefore,  a  new  expansion 
must  be  sought.  In  reference  (10)  the  initial  orders  of  mag¬ 
nitude  of  the  flpw  variables  in  region  4  were  established  as 
well  as  the  thickness  of  this  region,  which  is  OC£«  /.  To 
first  Order  when  the  N-3  limit  is  applied,  the  constant  stream¬ 
line  velocity  result  of  the  region  2  analysis  must  be  replaced 
by  an  expression  for  velocity  which  accounts  for  pressure  gra¬ 
dients  along  streamlines  which  are  found  to  be  the  same  for  all 
streamlines  in  this  layer.  After  determining  expressions  for 
the  flow  variables  in  region  4  matching  with  region  2  is  then 
demonstrated  and  a  composite  expansion  is  then  formed  which 
is  valid  throughout  regions  2  and  4.  For  general  shock  shapes 
these  expressions  must  be  numerically  integrated  to  yield  va¬ 
lues  for  the  flow  variables  of  interest.  Having  established 
results  for  distances  of  0(1)  away  from  the  axis,  we  now  wish 
to  determine  expansions  valid  near  the  axis  in  the  N-B  limit 
for  use  in  the  initial  radii  of  curvatures  of  the  shocks 
scheme . 

In  this  study,  unlike  in  past  investigations,10'  21'  22  atten¬ 
tion  is  focused  on  this  near  axis  region  because  we  are  seek¬ 
ing  accurate  and  simple  analytic  expressions  for  the  layer 
thickness  and  contact  surface  pressure  distribution  as  func¬ 
tions  of  shock  radius  of  curvature  to  be  used  in  the  initial¬ 
izing  scheme.  It  will  be  shown  that  the  near  axis  region  1 
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will  be  ofO(€«)in  thickness  andO(£4)in  length.  This  region  is 
characterized  by  constant  density  to  first  order.  As  in  the 
region  away  from  the  axis,  the  first  order  velocity  along 
streamlines  is  constant.  Consequently,  as  in  the  away  from 
axis  case  in  region  3  near  the  contact  surface  we  must  mod- 

•  t 

ify  "the  expansions  to  include  the  effect  of  non-constancy  of 
velocity  along  streamlines.  The  dimensions  of  this  region 
are  in  thickness  and  oW)  in  length.  The  variables 

found  in  this  region  are  demonstrated  to  match  with  those  of 
region  1.  A  composite  expansion  for  region  1  and  region  3 
is  then  formulated.  Finally,  matching  between  1,  3  and  2,  4 
is  demonstrated  for  the  specific  case  of  a  spherical  exter¬ 
nal  shock. 

For  the  internal  layer  analysis,  more  details  will  generally 
be  included  than  for  the  external  flow  since  this  work  ap¬ 
pears  here  for  the  first  time.  For  region  6,  which  is  Ofl)  in 
length,  the  N-B  limit  leads  to  the  results  that  the ( layer  thick¬ 
ness  and  streamline  velocity  to  first  order,  for  non-trivial 
results,  must  be  of  0(€t 'J2)  .  Also,  the  internal  shock  is 
spherical  in  form  and  the  first  order  pressure  and  density  are 
constant  across  the  layer  'with  an  error  0(€  «V*0  These  re¬ 
sults  will  lead  to  interesting  conclusions,  which  will  be  dis- 
cussed  at  the  end  of  the  analysis.  The  solution  found  in  6 
unlike  the  external  flow  case  is  valid  throughout  the  shock 
layer  for  distances  of  O(l)  .  For  general  shock  shapes,  these 
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equations  must  be  numerically  integrated  away  from  the  axis. 

As  in  the  external  flow  case,  a  region  of  0(£i J*)  j_n  thick¬ 
ness  and  of  o  in  length  is  required  near  the  axis  to 

be  used  in  the  initialising  scheme.  As  for  the  external  la¬ 
yer,  the  near  axis  region  5  is  found  to  have  constant  density 
to  first  order.  A  solution  technique  is  then  applied  which 
makes  use  of  this  characteristic  to  determine  analytic  expres¬ 
sions  for  internal  layer  thickness  and  contact  surface  pres¬ 
sure  distribution  near  the  axis  to  be  used  in  the  initializ¬ 
ing  scheme.  Having  found  solutions  for  the  variables  in  re¬ 
gion  5 i  we  then  match  these  to  the  variables  obtained  in  re¬ 
gion  6  for  the  specific  case  of  a  spherical  shock. 

The  detailed  analyses  for  each  of  the  six  regions  will  now  be 
given. 
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External  Layer  -  Region  2  Analvs-i  a 

Following  previous  investigations,21  we  obtain  the  orders  of 
magnitude  of  the  variables  inside  the  layer  from  the  values 
obtained  from  the  shock  relations.  These  relations  are: 


P* 

u 


A 

V 


P«e  =  J3«oe  ^«o«  ( 1 

U  ©*  COS  O' 

U»N  A 

toe  £  S MO  or 


-  t) 


a 

sin  c r 


/  =  l/i  -  Vi  ^ 


1 

y«+i 


Or**i>n«£si«*£ 


) 


(1.1) 

(1.2) 

(1.3) 

(1.4) 


Non  dimensionalizing  the  pressure  by  j^U^,  velocity  by 
Uooe.  and  density  by  ^  ,  we  obtain  from  (1.1-1. 4) 


p  *  y^eK «£  +  ( 1  -  Z)  si*)zc r 

A.  A 

U  ~  05 s  cr 

f')  A  A 

V  ~  e  sm  cr 

f  «  V? 

Taking  the  NB  limit,  we  find  from  (2. 1-2.4) 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


p  »  sin  cr 
U  s  Coscr 
0  &  e.esm.a' 

f  -  V€e 


(3.1) 

(3.2) 

(3.3) 

(3.4) 
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Using  the  orders  of  magnitude  established  by  (3. 1-3.4),  we 
can  now  write  the  following  expansions: 


U  * 
V  - 

f>  - 

r  - 


tlo  ■*■  £e  Ui  +  ■  • 
£eVo  +•  £e*  V*.  +  - 

P©  +  £«  pi  +  • 1 
jV««  +  ^a.  +  • 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


In  order  that  the  flow  variables  described  previously  lead 
to  a  meaningful  description  of  the  flow,  we  must  stretch  the 

normal  coordinate  y.  From  the  continuity  equation  we  find 
for  y 


“as  Yrx 


^u^rrr 


fixe  Uopf'rr — .  -v.  r€.e 
x  £^*-LUe  'X-TT 

ce 


Since  r  -o(RSJl)  which  is  set  as  the  unit  length  scale,  then 
y  — 0(<£e).  Therefore,  to  obtain  non— trivial  results,  we  must 
stretch  y  byO(€«).  The  expansions  of  the  variables  and  the 
stretched  coordinate  y  Should  now  be  applied  to  the  govern¬ 
ing  equations 


=  o 


H  . 


U*&U  -t-  H  v"^L  —  KUV  4-  J-  'b 


^X. 


J-  QJ3  - 

-P  Bk  ~ 


u^v  +  Hv'by  —  Ku1  h-  -H 


'bx 


j3 


X  -fc 
o 
o 


u.~ap/py  +  v  H  ^  P/f*  _  0 


(5.1) 

(5.2) 

(5.3) 

(5.4) 
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Before  substituting  in  the  expansions  of  the  variables,  it 
will  be  advantageous  for  final  integration  of  the  equations 
to  change  the  above  system  to  one  employing  the  stream  func¬ 
tion  as  an  independent  variable. 


Define  ' 

~av  _  u 

■Sx  -  H 

rpv 

(6.1) 

^  =  " 

rjaut 

(6.2) 

The  above  : 

system  under 

the  transformation 

=  ^ 

■&X 

■a  _ 

:  'ciV  d 

then  becomes 

H  r  jbfe.  = 

W  bx 

•  K/U- 

(7.1) 

Cu*  +v 

*)/z  ■+•  h 

■* 

’■+  V**)/*.  + 

h* 

(7.2) 

=  p*/p* 

(7.3) 

=  Hsi. 

(7.4) 

Bx 

u. 

The  distance  across  the 

layer 

in  terms  of  ¥ 

is 

obtained  from 

c)^  -  — 

5$ 

_J _ 

rjsa 

■ 

(7.5) 
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In  this  system  we  then  must  add 


r  =  rc  -4  £erx  +  •  •  •  (8.1) 

H  -  1  ~  t (8.2) 

Expanded  boundary  conditions  at  the  shock  are: 


/s 

p  = 

£e/M4€.  +  Sun1?  (9-1) 

A.  _ 

u  - 

COS  ^ 

(9.2) 

V  * 

Cl  +  Cot  ^CT  ^  Sluff- 

(9.3) 

rwi<* 

/V 

f 

1/  C  €•  (  1  +  )) 

«re 

(9.4) 

This 

then  gives  the  result  for  first  and  second 

order  boun- 

dary  conditions: 

AN. 

Am 

P»  ’ 

!  sm'ff 

(10.1) 

Uo 

=  cos 

(10.2) 

A 

Vfe  - 

:  (l  4*  coi^r/M^eeisiKi^o1 

(10.3) 

A. 

(10.4) 

Pi  5 

-  VM«tl£t  -  Ci  +  cot*<5/Moot€^  )sih*f 

(10.5) 

Qi 

=  O 

(10.6) 

-A. 

Vi 

-  O 

(10.7) 

A 

f1 

-  o 

(10.8) 

Substituting  (4. 1-4.4)  into  (5. 1-5. 4)  and  then  collecting 
terms  of  like  order  in  €e  we  obtain 
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First  Order  Problem 


T 


dy 


(11.1) 


fo  -  p®  po*  /  po* 

Uo  —  Uoj|t 

v-  *  3  Vs 

»•  ■  r  -is— 

U  Jvj,  /^Uo  To 


(11.2) 


(11.3) 


(11.4) 


(11.5) 


Second  Order  Problem 


r 

Pi  =  fi-  rtSsUi  +J-5&  +  +  a.«s<r)]<l'p  (12.1) 

r  Js<  r«  r»ax  r‘  (12<2) 

fl  =l\pi*/F>3.  +(PVp.  -  P»*/fc*)jV  h/posr  ‘ 

111  =  ~  ■£»*_  I M  P>/Pos»  , 

UoPo*  (12.4) 

Vl  *  Uo^-  +  KaJf»v«  +  v«Cu»/u.P  (12  5) 

1%  *  (WP-  +  Ul/u»  +  ri/r->) 

Equations  (ll.l)to  (12.5)  describe  all  the  properties  of  in¬ 
terest  near  the  shock,  with  an  error  o(e.Vaj. 


(12.3) 

(12.4) 


(12.5) 


External  Layer  -  Region  4  Analysis 

Bush  10  ha?  determined  from  these  results  that  equations 


(11.1)  to  (12.5)  are  rot  valid  at  the  contact  surface  due  to 
the  change  in  order  of  magnitude  of  streamline  velocity  and 
hence  another  expansion  is  required.  Following  his  results, 
we  derive  the  system  of  equations'  valid  in  the  contact  sur¬ 
face  region  with  an  error  0(e*j. 


1 
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The  expansions  for  the  variables  in  this  region  are: 


IL  —  €e  Uo  +  6ft  Ui  t  •* 

(13.1) 

V  ®  Se*  \£  +  6^  Vl  4-  •  •  •  • 

(13.2) 

P  ~  p©  +  pi 

(13.3) 

+  ft  +  *  *  ' 

(13.4) 

l"  -  re  4  €eVa  rx  4*  •  •  * 

(13.5) 

H  -  1  +  KLo  Jo  ecVa  + - 

(13.6) 

7C  =  5c 

(13.7) 

3-  " 

(13. B) 

These  will  be  substituted  into  the  governing  equations  (5.1) 
to  (5.4).  The  resulting  systems  of  equations  ares 
First  order  system  of  equations 


Second  order  system  of  equations 


(H.l) 

(14.2) 

(14.3) 

(14.4) 


(15.1) 

(15.2) 

(15.3) 
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We  can  extract  from  (14.1)  to  (15.4) 


First  order  results 


P°  *  P“°°  (16.1) 

■t?0x  =  UoJ  -  zpiyjS,  In  Po/pi*  (i6.2) 

fo  -  fT.jWp^,  (16.3) 

Second  order  results* 

Pi  *  ^  (17.1) 

X.U0U1  =  ZMcUx*  ~  C  ^  “  (17.2) 


(17.3) 

+ (fe*  ~  W 

Po*  r*»  r°* ' 

These  body  layer  results  must  be  matched  to  the  external  la¬ 
yer  region  2  results  in  an  intermediate  region  where  both  ex¬ 
pansions  are  valid. 


*  is  not  considered  since  it  does  not  affect  the  calcu¬ 
lation  of  P ,ML.  to  error  greater  than  ofCe^and  is  itself  very 
small  of  OC€«*-); 
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—  -w .i.. if.. 


We  now  demonstrate  matching  for  the  variables  and  U.  . 
Matching  occurs  within  the  intermediate  region  defined  by  7\  . 
In  this  region  . 

O  *5  — ►  oo 

In  terms  of  the  variables  <$>%  and  9  j^withf,  ^  of  0(l)  we  have 
from  the  expression  for  layer  thickness,  (7.5), 

y*  -  ou\  <p*  =  o(^1) 

The  order  of  ?{  must  then  fall  in  the  range 

o(€^)<  ^  <  o  (i) 

2  3 

Following  Cole  ,  matching  will  be  achieved  if  for  each  power 
of  we  have 

\\yy\  Q  ?  "  f  ^  O 

€e-vO 

where  f  is  the  variable  of  interest  to  be  matched.  Taking 
the  matching  of  the  pressure  first 

Wjj  Cpc?)  -  pC$.,$>*-*0)) 
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Therefore,  the  pressure  in  the  body  layer  is  equal  to  the 
pressure  at  the  inner  boundary  of  the  external  layer.  A  com¬ 
posite  expansion  valid  in  both  regions  can  be  determined  by 
subtracting  the  common  term  to  both  expansions  which  is  the 
pressure  at  the  contact  surface.  The  correct  pressure  dis¬ 
tribution  will  then  be  that  derived  from  the  external  flow. 
For  matching  of  the  density,  we  must  have 

From  (16.3)  and  (17.3)  we  have 


Since  P«  and  pi  are  just  functions  of  thenp(®)  ;  there¬ 
fore 


U  (  pcf )  -  f(Q,o))  $0 


The  density  within  the  body  layer  is  then  equal  to  that  at 
the  inner  boundary  of  the  external  layer.  As  in  the  case  of 
the  pressure  distribution,  the  uniformly  valid  density  distri¬ 
bution  is  simply  given  by  that  derived  for  the  external  layer. 
The  velocity  now  has  to  be  matched.  Once  again,  as  for  p  and 
P  ,  we  must  find 
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Iirv* 

€e-»0  ) 


For  the  body  layer,  we  have 


Uo*  = 


Pb*  “  ft 


o* 


=  1 


U*.*  =  O 


pi*  =  '1  - 


~  _ 
f  4*  *  ° 


Under  the  above  conditions  U  from  (16.2)  and  (17.2)  becomes 

u  =  €.**■(§*  -  2.ln'poC9)Vi  + 

c-  vmK  m -vm-fe.  -  A)  - 

(e*  - 


The  near  shock  expression  for  the  velocity  is  from  (11.3)  and 
(12.3) 


u  -  0*  -  fog,  In  pb/jp^ 

O*  jS* 


For  matching  in  the  intermediate  region  ^  we  transform  the 
independent  variables  such  that 

©*  -  CO  ^  O 

e*  =  e*7p\  =*.  o  <$>*  =  <p»v?i  *  o  ea  *  o 

34 


The  matching  condition  then  becomes 


linn  (  Q*ti71  —  lo  po  — .  C  In  po)  6a  +  •  •  • 

6e*ox  v  e.1'*  e„n7i 

+  £e.  (  ~  PV  fi.  Ml4^e  ~  *•  ~  (.VflajSe  ~  0  ^ 10  P°) 

(&2L-& ^  •  •  ) 

e*Vl-  1^71  2.64,71* 

-  0*7)7]  ■+  €»  In  pe(fco)  )  =>  O 

©*7,?l 

The  above  can  be  seen  to  match  to  error  O ((■  e^*)  .  A  compo- 
site  expansion  valid  in  both  regions  is  then 


Ce**  -  z«£elnpoC^,o)) - 
€e  (  P°Vja»^n  P°/po»  ~  'n  P*Cf'°)) 

(©**  -  26eln  p0C9/0))1/a- 


(18.1) 


The  layer  thickness  can  now  be  calculated  from 


Once  the  radius  of  curvature  and  position  of  the  shock  are 
known,  then  all  the  layer  properties  may  be  calculated. 


Now  that  we  have  the  results  for  the  external  flow  away  from 
the  axis,  we  will  now  determine  if  they  apply  in  the  region 

near  the  axis.  'From  the  external  flow,  the  velocity  parallel 
to  the  shock  is  given  by 

U  =  6*  -  -po  c  9,0) 

As  we  approach  the  axis  both  $  and  <J>*  ,  as  well  as  6*  be¬ 
come  smaller  than  0(1)  .  For  this  situation,  the  expansions 

based  on  the  above  quantities  being  of  order  one  are  then  no 
longer  valid  and  new  expansions  must  be  obtained.  It  is  of 
importance  to  develop  an  accurate  solution  in  this  near  axis 
region  as  the  final  matching  of  the  external  stream  layer  flow 
and  jet  layer  flow  depend  on  knowing  the  initial  locations 

and  radii  of  curvature  of  the  shock  surfaces  near  the  axis  of 
symmetry.  V 


External  Laver  -  Region  1  Analysis 

To  establish  these  initial  values  for  the  external  flow,  it 
is  convenient  to  determine  the  shock  layer  thickness  and  con¬ 
tact  surface  pressure  distribution  in  an  analytic  form.  This 


leads  us  to  describe  the  external  shock  as  an  expansion  in 

even  powers  of  <J>  for  the  region  near  the  axis.  The  expan¬ 
sion  is 

R,  -  1  4-  9  *4  be  (p  ^  4  *  • 
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It  is  sufficient  to  retain  only  the  first  two  of  the  above 
terras  for  accuracy  consistent  with  the  number  of  terras  in  the 
dependent  variable  expansions.  The  expansions  for  the  de¬ 
pendent  variables  which  give  the  relevant  physical  behavior 
in  the  near  axis  case  and  are  consistent  with  the  shock  re¬ 
lations  are 

p  ~  p©/  ■+■  £e.  pi  +  •  ' 

f  '  po'/e*  4  ft  4  ’  *  *  ’ 

U  -  U0'  +  £eV*  Ul  +  *•’ 

V  -  Vo'  +  Vi'  +■  •  *  • 

Inserting  tho  above  in  equations  (7.1)  -  (7.5),  we  obtain  for 
the  first  order  results 


For  the  second  order,  we  obtain 


(19.1) 

(19.2) 

(19.3) 

(19.4) 

(19.5) 


(20.1) 

(20.2) 


(20.3) 

(20.4) 


ZUo  Ui  +  Vo'  *  =  ^u4  ui*  +  Vo*  -  X. 


The  boundary  conditions  obtained  from  (1.1)  -  (1.4)  are 


-  1  Uo*  =  (?*  (A  “  Vo*  =■  1 

pj.4  *  -  (l  +  ~  VmwSc.) 

fti.  -  -  Ui*  -  -  9i*C l-Xae)5/fe 

00c 

For  R-»l  +  a*q)  (19.1  to  19.5)  and  (20.1  to  20.4)  give  for 
first  order 


fb'  =  fo'  =1 

Uo‘  =  Cl  -  2-a«  )  9*' 

Vo'  =  (<P*'/§')*' 

!t‘  -o^55U  ■ 

For  second  order,  we  have 


(21.1) 

(21.2) 

(21.3) 

(21.4) 


pi  •=  -(l  +  $'*(l-2fl«)*'  -  l/M«6e)  — 

((i-2a  ef  9^r)/3  -  (i-^-rVa) 
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External  Laver  -  Region  3  Analysis 


As  in  the  case  away  from  the  axis,  we  must  develop  a  body 
layer  solution  near  the  contact  surface.  The  relevant  ex¬ 
pansions  for  the  variables  in  this  region  are 

p  *  +  £  e  pjJ'  +  .  .  • 

J°  “  f°"/ +  fl  + - 

LL  -*•  €e  -f 

V  -  de*"  Vo''  +•** 

X  -  6.x  *  =  6.  £ 

Substituting  the  above  in  the  equations  (7.1)  to  (7.5),  and 
collecting  terms  of  like  power  in  €«  »  we  obtain  for  first 
order  results 

^  *  O 

(23.1) 

PoVfc"  =  pA»/f>oi  (23.2) 

For  second  order  results 


ii 

o 

(24.1) 

Uo"  1  =  U.V-  4-  zCPil  - 

Pf') 

(24.2) 

fi"  *£*±tL  +r£*T-Bi\ 
Poi'  ^Po"  p.;> 

pJLe^' 

R>» 

(24.3) 

The  boundary  conditions  which  must  be  applied  are 

po*  •=  J^'  -l  Uo*  =  (l  ~  a«e)  9* 
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- - - — 


pi*  -  VM&e.  -t 

% 

Ths  final  results  which  must  be  matched  to  the  near  axis  outer 


region . case 

are  then 

n 

po 

=  po  lx) 

(25.1) 

_  1/ 

«  />."<*> 

(25.2) 

^  II 

p* 

=  PiCx) 

(25.3) 

_  II 

p' 

=  A"c*> 

(25.4) 

p±  =  JVC*) 

Uo"  =  CO.  -*uf  #*+  2-Cpi'i  -  pO)  (25-5) 


Following  the  same  procedure  as  in  the  away  from  the  axis 
case,  matching  requires  that 

peCX)  ■*  po'  =  J. 

po'Cx)  =  pa'  =  1 

p"  =  pk§,o)  =  -  Ci 

pi"  =  fi'C$,0)  =  -*>«)*  +  i; 

Th“  values  obtained  for  (26.3)  and  (26.1)  are  simply  those  de- 
rived  from  the  near  shock  region  expressions  evaluated  at  the 
contact  surface.  The  value  of  U©"is  then 

Uo"  =  (Cl-aa.)^;*-  +  £(-£+^$'*(l-2«e)l)(26.5) 


(26.1) 

(26.2) 

(26.3) 

(26.4) 
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Matching  to  outer  layer  results,  we  require  that 

\\*rt  (  U*  -  U")  $  O 

0 


:e+Ov  '  aa-2^0 

(- f ' 1  -  9n?)Zp j/3 - 

4(1  -  *W®  -  **  (( + 


Expanding  the  square  root 


€*-*0 


%ZZq>')/3  - 

-^(l-  9*  V/ ^,y)  -  £e  (&-*«.)& VeeVa  + 

1  -  *»«)*•  -  V^Kl-iaO^n  - 

♦.£&£*))*> 


**  From  boundary  conditions  derived  on  the  basis  of  mass  flow 
through  body  layer 


Tt  is  seen  that  matching  is  assured*  A  uniformly  valid  re— 
suit  for  velocity  in  both  regions  is  then 

1/2 

U  =  (  6^((l-3.a.)l^+ 

+  6ev*  2C-f 1 *a-*«o*+ct-a^‘(r - 

5>is/f)/3  -i(i  -  <W0) 

^  1  -  W.)‘^+  ^  £c ($$ rf\  27.1) 

Using  the  uniformly  valid  results,  we  calculate  the  layer 
thickness  from  (22.2)  (27*1)  and  (7*5),  which  gives 


+■  e*x 

(l-aae) 


+  «&-*«•>(&■«■&  +3fe)f'  (28-1) 


The  contact  surface  pressure  distribution  is  given  by  (21.1) 
and  (22.1) 


Pcs  ~  1  ~  -  (1-  ZQ.')Z±<p‘  €e  +  _£s -  (28.2) 

r  a  3 7  M«ie« 

Expressions  (28.1)  and  (28.2)  are  then  used  in  a  scheme,  out- 
lined  in  Appendix  A,  which  determines  the  external  shock  ra¬ 
dius  of  curvature  and  position  along  the  axis  of  symmetry, 
for  given  external  flow  and  jet  conditions.  These  values  are 
then  used  as  initial  conditions  for  determining  flow  proper¬ 
ties  away  from  the  axis. 
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The  solutions  found  for  regions  1  and  3  as  the  variable 
P  tends  to  infinity  must  match  the  solution  found  for  re- 
gions  2  and  U  as  vanishes.  If  matching  is  achieved,  then 
we  have  proof  that  the  expansions  for  the  region  near  the 
axis  was  the  correct  one. 


For  this  analysis,  a  composite  expansion  lor  the  near  axis 
and  away  from  axis  regions  is  not  sought.  Therefore  a  match¬ 
ing  of  these  regions  for  the  special  case  of  a  spherical  shock 
will  demonstrate  the  consistency  of  the  expansions  and  will  in¬ 
volve  a  much  smaller  calculational  effort  than  if  the  general 
shock  case  were  attempted. 


The  away  from  axis  results,  equations  (12.1),  (12.2),  (18. l), 

>s 

for  a  spherical  shock  and  small  are 

=  <p,*  -  ?*73  +  *e= (§^-9*739-  9*7 


p>  -  1  -  €«  - 

r  Z 


37  3©  I*-7 


tom 


3  12  ^s.  3  If  a#? 


3  W 


~ 

2 


j* =  1 


-  a.8  -  j-9*m  +  x 


+  -  9*1  +  (1  +  9**~ 

3  3  sgf  v  m»;«. 


-V)  +  ^l(  1  -  Tfe-jX1 
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»  v*U4  W-tt  £<m 


+*  l>nio»<£iaVC«UK'..  ■**> a-  W. 


MrifeMifiiiiiiHiiiiMiiiiM 


The  near  axis  results  for  a  spherical  shock,  Cl«-0  »  from 
equations  (27.1),  (21.1),  (22.1)  and  (22.2)  are  given  by 

U*  =  £e  9#'1  +  -z(-$‘z+  8'*- &'*/#') 

-  (i 

P  “  1  —  H/3  €e  9  +  £e9#/3^  "+  £e/M«|fe 

- 

J=  =  1  +  +  etOf-*1-?'4)  ~ 

9,t-<fr,9y 


Following  the  same  matching  procedure  as  was  outlined  for  re¬ 
gions  1  and  3  and  2  and  4,  we  define  a  variable  7^  which  des¬ 
cribes  an  intermediate  region  where  both  the  away  from  axis 
and  near  axis  expansions  are  valid.  For  this  case,Ofe^  ?^<0(l) 
and  therefore  for  li  matching 

5§n7l 

(fe'  - tg$)  * 


Matching  is  assured  for  U.^*  . 


t 


For  p  we  have 

Uo(i  -€«/*.  -^9„V/5  +  q>*%7lY3^n  - 

~  +■  £e/€eM»«  -  £<s<P*q/x.§t  + 

T}V3%’  -  €e.<p*n7l*/zf7)i  -  1  -h  €e/x.  + 

4^lV/3  -  -  6^4ni,  + 

«•$***  A  )*0 


Therefore,  matching  of  p  to  this  order  is  assured. 

For^p  we  have 

tlTo  ^  1  '  +  9*W/sfa  -  9*« (Vm.v*  -  l) 

+  €c  (- 1  +  7)l<^'(l-  1/M<f«€«  -12/12.)  +  1/2-  + 

4&9M*/*3  -  -  $>*yr/3?„  -g£n/z$J 

+  *%UV-39n  +  e=/2.  -  <^&y<pn4  -  1  - 

9*nflVM«*6e  -  (9*n  ~  <ph)  ?? 2  +  C fl*/3  “ 

-€eA(  1  -  %\/&)))*0 

Therefore!  matching  p  to  this  order  is  assured. 
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Internal  Layer  -  Region  6  Analysis 

We  start  the  analysis  from  the  shock  relations,  which  are: 


p  =  p»t  Mooi  (l  ~  ^ ) 

_  ■O'  y  A 

f3  -  /><*;  /  e 
u  *  U®i  cos  & 

A,  /V  I  |  A 

V  *  €  Uooi  sm/5 

£  =  *:-i  4 

sy*i 


(29.1) 

(29.2) 

(29.3) 

(29.4) 

(29.5) 


Unlike  the  external  uniform  stream  case  fa,  fai  are  not 
constants  but  functions  of  the  undisturbed  plume  flow.  Their 
functional  relationships  are  determined  by  the  assumed  source 
distribution 


We  now  expand  the  boundary  conditions  for ^O/V-kjo  such  that 
O(l)  .  This  will  then  indicate  the  orders 

of  magnitude  of  the  variables  within  the  layer  near  the  shock. 


It  is  hoped  that  the  orders  of  magnitude  found  for  the  vari¬ 
ables  will  be  valid  through  to  the  contact  surface.  If  this 
is  not  the  case,  then  the  solution  based  on  the  orders  of  mag¬ 
nitude  near  the  shock  will  break  down  and  a  second  expansion 
valid  near  the  contact  surface  will  be  required.  Expanding 
the  undisturbed  plume  parameters,  we  find  that 


When  applied  to  (30.1) , (30.2)  and  (30.3),  we  obtain 


(31.1) 

(31.2) 


(31.3) 


The  first  order  expressions,  (31.1),  (31.2),  and  (31.3),  are 
then  introduced  into  the  expanded  form  of  the  Rankine-Hugoniot 
relations,  (29.1),  (29.2),  (29.3),  (29.4)  and  (29.5),  which 
then  give  the  order  of  magnitude  results 


f  ~~  °(£o) P  O(t)  fOD^UitoC 
U  ~  0(.©)  Uooc  V  0(^ce)  Llooc 


where  ©  is  the  angle  the  free  stream  makes  with  the  direction 


normal  to  the  shock. 


At  the  axis  of  symmetry,  ft'o  is  the  density  ratio  across  the 
shock  and  is  given  by  6J©  =  ¥.t  ~J=-  -+-  jjT  .  Up  to 

ifu4-  1.  »"«Ol6 

this  point,  the  magnitude  of  the  velocity  u  parallel  to 
the  shock  has  not  been  determined  in  terms  of  6i0  ;  however, 
its  lower  bound  can  be  determined  from  the  condition  that  the 
shock  layer  is  of  zero  thickness  in  the  limit  as  and, 

therefore,  the  flow  must  be  turned  90  degrees  upon  crossing 
the  shock.  The  deflection  angle  is  given  by 

=  i  s\nx&  -  l)coj^ - 

*5^-  (P)o£s«"»i5  “  1) 

©  for  near  normal  shock  © <<  1  Moot' ^  1- 

"ban  §  —  cos*Q  ban  0  -  © _  ^  00 

(***  -  cos*ej 

This  then  yields 

lnrv^  Q _  ^  co 

6U  +  0  +  ©*) 

We  can  conclude  that  ©  ^  0(€t'o)  and  is  obviously  less  than 
0(1)  .  This  result  can  be  written  as 

U  =■  0(€co  )  l^hn>0 

The  dependent  variables  can  now  be  expanded  in  terms  of  ©l'o  » 


•ban  S 


where  -  JL  ~ 
7L 
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JLi 


following  the  first  order  magnitudes  found  from  the  shock  re 
lations.  The  non-dimensionalized  expansions  will  then  be: 


p/jf^Uool- 

~  P>0  + 

+■  Pa  +  ’ 

f/fooio  = 

jPo/  Slo 

i  +■  + 

U/  Ucoi  ~ 

U0 

+  +  6  u 

V/  UoOt’  = 

e;„  Vo 

It  is  assumed  that  X  will  be  of  order  one,  from  geometry, and 
will  depend  on  to  some  power  n*  The  unknown  pow¬ 

ers  m  and  n  will  be  determined  by  requiring  that  the  govern¬ 
ing  equations,  when  written  in  terms  of  £t0  ,  yield  a  physic¬ 
ally  reasonable  non-trivial  system  of  equations.  The  govern¬ 
ing  equations  are: 


continuity 
X  momentum 


momentum 


energy 


where  K/  is  the  radii  of  curvature  of  the  shock  surface  and 


H  “  is  che  ratio  of  the  radius  of  curvature  of  a  constant 

surface,  with  respect  to  the  shock  curvature  at  a  given 


value  of  X. 


Substituting  the  expansions  into  the  governing  equations,  we 
obtain  the  following  results  relating  the  unknowns  m  and  n*: 


■te*  =  oter-1) 

t:  o(ecf  )  -  ©(sj."1) 
£i<,m  =•  o(€»l‘") 


X  momentum 
g  momentum 
Continuity- 


Applying  the  source  condition  0(l)we  find  that  m  =  i. 

oX 

It  then  follows  that  0  =  a  from  the  other  two  relationships. 
Since  6*^-^  is  now  found  to  be  proportional  to  €  , 

then  in  the  limit  as  €t*o  0  ,  must  equal  7T/£  .  The 
shock  surface  is  then  a  spherical  shell  which  is  centered  at 
the  jet  exit.  The  magnitude  of  the  perturbation  from  this 
spherical  shell  will  then  be  derived  b]  the  use  of  the  known 
value  of  ©  -  0(6to X) ,  The  radius  can  be  written  as  R=l-«-GCx) 
From  this,  the  value  of  ©  can  be  calculated  R'(x)  =  G'oO  =  6^0 
*  The  magnitude  of  Stx)  is  then  o(e  <.'•*)  which  then  leads  to 
the  result  F«X>  =  1  +  FOO  6t'.Va-  .  This  will  be  the 

form  used  for  R(x>  .  With  the  values  of  m  and  n 'established 
the  expansions  are 


p  ’ 

Po  +  pl 

Guf*  +  *  *  *  * 

f  ’ 

=  J*  o/€co 

+  fn/ec 0V*  + 

U  ■ 

-  ^.V2  u 

i  U  x  "h  €  Mid*  • 

V 

-  See  Vo 

+  Vi  +  €t<0*  vx  + 

* 

_  X  ,  Vx  u 

-  6*  y 

X  -  x 

The  flow  equations  become  considerably  simplified  for  inte¬ 
gration  across  the  layer  if  they  are  written  in  terms  of  the 
Mises  variables  X,  V  . 


The  distance  downstream  is  X  and  V  is  the  streamf unction. 
The  flow  equations  in  this  coordinate  system  are 


Hr^.  =.  K.U  d- 

^i±Xl  +  -£il  -  COfO 

2.  *rx? 


?/f*  *  blv) 


H  -  l  - 


(32.1) a^  momen¬ 
tum  equation 

(32.2)  inte¬ 
grated  X  mo¬ 
mentum 

(32.3)  energy 
equation 

(32.4)  stream- 
function  defi¬ 
nitions  which 
satisfies  con¬ 
tinuity  direct¬ 
ly 


(32.5) 


£fter  substitution  of  these  expansions  and  introducing  an  ex¬ 
pansion  for  ^  and  K  in  terms  of  X  and  4*  and  collecting 
terms  of  like  order,  we  then  develop  the  following  system  of 
equations. 

First  Order  Expressions 
Momentum  equation  across  layer 

-  o  (33.1) 
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Integrated 


momentum  equation 


Uo*  -  Uoi  -  In  Pft.  (33.2) 

fo*A  pb* 

Integrated  energy  equation 

po  ~  fo  (33.3) 

Definitions  of  streamf unction  equations 

i -  (33.4) 

c)  rQ 

Definition  of  streamfunction 

-  ¥s.  (33.5) 

Uo 

For  the  second  order  results,  we  obtain  for  the  momentum 
equation  across  layer 

"Spy  -  K.o  Uo  (34.1) 

ro 

Integrated  momentum  equation 

2UoUi  -  2.Uo*Ut*  -  3.po»  Py  (34.2) 

?o 
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Integrated  energy  equation 


pi  =  f  i 


(34.3) 


Definition  of  streamf unction  equations 


1 

^oUofb 


+ 


Ux 

Uo 


+ 


(34.4) 


Definition  of  streamfunction 


^4i  -  Kc4.oVft.  -+-  W  -  Va  Ux  (34.5) 

°  Uo  Uo  Uea. 


For  the  third  order  results,  we  obtain  momentum  equation 
across  layer 


-  KqM3x  +  ■+"  UflJs^rokcifo-rL)  n,  r 

r0  >o  fcx  «7,«-  ^ 


Integrated  momentum  equation 


Ui*  +  ZUzMo  +  V0*  =  2-Ua*U0*  +  Mi*4  +  (35.2) 

L£&t  lh  £s. 

J~k*t  R* 

"  %)  _  *-ft  llh  fJ 


Vo*  -f  3-Pt >»/>*»  In  Po  -  jLfiZff  lh  js.  — 
/*** 

\SL 
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4u^.i — J.f  -l.'— 


Integrated  energy  equation 


Definitions  of  streamf unction 


- - 1 -  (  -  (gi  +  U» 

jo0uo  n>  \  u* 

UvJDt^  +  (£±  +  3.  +  tu 
u.6ro  y  \go  ni  Uo 


r) 


-+■  fia. *jl  +  ga  Ut 

J^o  *"0  J>o  U-o 

(35.4) 


We  must  apply  the  expanded  shock  boundary  conditions  to  the 
outer  edge  of  the  layer.  In  the  limit  as  6 1'0  “►  O  ,  fico^-^oo 
such  that  tco  Meo^0  OC1)  we  obtain 
First  order  results 


p° 

.A. 

-  cosn m 

X©oo 

(  36.1  ) 

J5* 

"  COS*  TTCP- 
2l©co 

(  36.2) 

A. 

Uo 

=  F'cf) 

(36.3) 

Ao. 

Vo 

-  l 

(36.4) 

a* 

-  0 

(36.5) 

Second  order  results 


P* 


-  ^Fc«"; 
-ZF  cos" 


>00 


(37.1) 

(37.2) 


U»  =  -  F'F 


Third  .order  results 


(37.3) 

(37.4) 


P1= 

as"^(3Fl-l 

-f/2*+  \  ) 

(33.1) 

_p*  = 

cos"M.  (.3F*  ■+ 

2 Bto 

J — — —  ^  in 

(&~OM 

-Tu  )  X 

00('o 

-  F/a) 

(33.2) 

Ui  = 

l  FF'  -  F'Vfe 

(33.3) 

>»x 

=  o 

(33.4) 

o 


Before  we  integrate  the  equations,  we  introduce  another 
transformation  of  the  independent  variable  V  .  This  vari¬ 
able  describes  the  streamfunction  in  terms  of  the  position 
at  which  the  streamline  has  crossed  the  shock.  It  is  de¬ 
fined  from  the  source  conditions  as: 


1^9*  sm9* 

2.6® 


(39.1) 


A  point  within  the  layer  will  then  be  described  by  the  co- 
ordinates  (  J>  >  <¥ *  )  .  Substituting  (36.1  to  39.1)  in  the 

equations  (33.-  to  35.4)  and  integrating  across  the  layer 
we  obtain 


First  order  system 


(40.1) 

(40.2) 

(40.3) 

(40.4) 

(40.5) 


Second  order  system 


<P 

pi  -  -2.Fco5n^o+  ±k  fi 


Ut  =  ~  aP^P*) 

f*  = 

.  /'T  *  * 


(41.2) 

(41.3) 


=  & +  k)  (w-° 

Vl  =  +  V^)  <41,-5> 


Third  order  system 


f>Z  -  foCpi/po  -  pi.fl/po4  +  (pa./p0‘)*'  - 
(pa.*/pe*  -  f **/fo»  +  ^Jrln^.))  (42.2) 

“  V^Uo  +■  Ui^  -h  Vo**  —  Uil  —  VoZ  4* 

(^F^/jPo»)fa^UjV^b»  -  In  p</p*-(42.3) 

*f*Vfa*OVpo  -P-V^)  - 


$ 


=  Vr.p0  f  «*"l||ksmg>*  Jg>*  (-(^  +  + 

Jg*  rrr - - 


Uo 


&-  +  Jgi^k  4-  RiUy  4-  UihA  4-  /£&.  +  Ui  + 

s,_A°  >U°  U°rJ  ^  ^  (42.4) 

•)  ) 


Hi 

r© 


Examining  the  above  integrals,  we  note  that  the  calculated 
variables  remain  finite  and  non-zero  from  the  vicinity  of 
the  shock  through  to  the  contact  surface.  This  indicates 
that  the  solution  is  uniformly  valid.  V/e  now  must  exam¬ 
ine  the  behavior  of  the  solution  as  we  approach  the  axis 
of  symmetry.  This  is  important  since  in  order  to  evaluate 
the  matching  of  the  shock  layers  we  must  determine  the  ini¬ 
tial  radius  of  curvature  of  the  shock  surface,  as  well  as 
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its  position  in  space.  This  requires  that  we  know  both  la¬ 
yer  thickness  and  contact  surface  pressure  in  an  analytic 
form  so  that  the  matching  of  contact  surface  pressure  and 
position  can  be  determined.  For  the  case  when  Rex)  -  i‘0 
the  variables  can  be  expanded  in  even  powers  of  near  the 
axis,  and  then  integrated  across  the  layer.  These  results 
will  then  be  used  to  illustrate  the  change  in  order  of  mag¬ 
nitude  of  some  of  the  dependent  variables  near  the  axis  of 
symmetry. 

From  the  source  conditions,  we  can  write  for  small  ^ 

P*  -  5  1 

which  then  leads  to  the  first  order  results 

P=  -  j>0  =  x  - 

=C^a(r-9w‘))ya 

V  =  UWC1  -M)a))'/a 

v.  =  f 

The  second  order  results  then  become 

P1  =  $&.($*  - 
f  «•  -■#•$.(? 1  - 

Ul  =  */s§  (f  *  -  9**) 


It  can  be  seen  that  as  we  approach  the  axis  -*■ 

and 

U  =■  €  K.  -s  0(0 

V  -  £U  Vo  -  O  (£,'.) 

p  -  po  +  px  =  a.  +■  o(els) 

j°  ~  fo/eco  +  P  i/e  =  oCVe.O  4  0(1) 

x  -  0(0*)  =  oCeJ1) 

The  change  in  orders  of  magnitude  of  the  variables  then  ne¬ 
cessitates  that  a  new  expansion  procedure  be  applied.  To 
first  order  p,j>  and  V  do  not  change;  however,  U  and  .X 
are  altered  and  the  new  expansions  must  reflect  this. 


Internal  Layer  -  Region  5  Analysis 

The  shock  position  in  this  region  is  described  by  an  expan- 

A 

sion  in  terras  of  even  powers  of  ^  away  from  the  axis;i.e., 


R  =  i  +  a,  6t*o9*+  • 


Based  once  again  on  the  shock  relations,  we  can  expand  the 
dependent  variables  in  the  form 

U  -  £toUo  +  +  '  '  ' 

V  =  +  €io3  W+  *  •  * 

p  =  po  +  €topJ.  +  *  *  * 

f  =  P°/e-  +  J*  +  ‘  ‘ 

—  V*.  —  -  V2l 

X  “  X  £u>  ^  ~  ^  0 


Substituting  the  above  in  the  governing  equations,  we  obtain 
the  first  order  system 


=  O 

=  o 


These  are  of  the  same  form  as  the  constant  density  flow  equa¬ 
tions. 
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uiSA-ito, ***3*0*1  ^ 


V/e  therefore  apply  a  constant  density  flow  solution  method 
to  this  first  order  system.  A  streamfunction  is  defined 
which  will  satisfy  the  continuity  equation  identically.  This 


is 

=  -FoUo  ^fc!  =.  r.Vo 

If  S-'  is  known  throughout  the  layer,  then  Uo  and  Vo  can 
be  calculated.  The  equation  for  S'*  to  this  order  of  approxi¬ 
mation  is 


VJhere  is  the  vorticity  within  the  layer.  For  an  incom¬ 
pressible  flow,  the  circulation  must  remain  constant  through¬ 
out  the  flow  region  and  this  leads  to  the  constancy  of  the 
product  of  vorticity  and  distance  away  from  the  axis.  Not¬ 
ing  this,  we  can  then  develop  the  vorticity  inside  the  layer 
in  terms  of  its  value  at  the  shock  and  its  radial  position 
within  the  flow.  The  result  is 

The  equation  that  must  be  solved  is  then 

%(f)  ^  -  -("fey*1'*-)?”'1' 
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Substituting  in  a  st re amf unction  of  the  form 


+  b^j. 4-  Ej£a  +  P  + 


(44.1) 


we  find  from  (43.1)  that 


A— •"<«) 


From  the  definition  of  the  st re amf unction  we  have 


VOL  = 

iBL  = 


=■  rov0 


=  -  To  Ut 


(45.1) 

(45.2) 


From  (45.1)  we  find 


-(n^*-  +  HqAq)^oa'  ZCQ  +  J.69^.  =  To  V0 

(46. 1) 

_  J_ 

At  —  O  which  is  the  shock  surface,  we  have  Vo  =  L 
which  then  gives  from  (46. 1)  C  “ 

From  (45.2)  we  find  —  ^  +  = 

at  the  shock  =  O  U0  *  and  from  (47.1) 

we  obtain 


6  =  0 


G  —  —  <x  i 


These  results  then  give  us  the  first  order  velocity  compo¬ 
nents  in  this  region. 

vo  =  1  -  -  (K§wT+‘<a‘)i’*'  (43.1) 

u.  =  ?  zai)  (48-2) 

and  po  -  po  - 

The  normal  coordinate  ^  to  first  order  as  a  function  of 
and  ^  can  now  be  calculated  from 

=.  JS-  (49.1) 

Ue 

Substituting  in  Vo  and  Uo  in  (49.1)  and  integrating,  we 
obtain  after  applying  the  shock  boundary  conditions 

1  -  +  -(§*/§) 

or 

o 

The  pressure  field  can  now  be  calculated  from  the  known  ve¬ 
locity  field  given  by  (48.1)  and  (48.2).  We  must  integrate 
the  strearawise  momentum  equation  to  obtain  the  variation  of 
pressure  along  the  streamlines  with  velocity.  This  result  is 

Uo  +  Vo*  =  Uo2  *  V„:  -  z(pi  -  P*»)  (51.1) 


s-CKifc'f +  40 


- ^fco.i) 
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From  the  shock  relations 


pi*  ~  C  1 


*ioc*-o 


Applying  this  to  (51.1),  we  then  obt; 


-  -  Vz  +  i/M«5,€f0  <p/2  ((h&)x+t/q^o  +  2.qj)' 
~(i  -  (52.1) 


From  the  conservation  of  entropy  along  streamlines  relation¬ 
ship,  we  also  have 


ft  “  pi  -  fi*  +  ft.*  (53.1) 

Making  use  of  the  value  of  pi-p^and  noting  from  the  shock 

relations  that  -  (  ft  **)  we  can  then  write 

for  (53.1) 

*  ***■*&*  -  C1  -^>-("^^0^)/^ 

-  +  *»;)%•  +  Va.  ~ 

(  ***)  9*a 

The  complete  expansions  for  the  dependent  variables  in  the 
near  axis  region  can  now  be  written 

u.  =  6 to  Uo  +  0(eie?) 

v  =  6^Vo  +  Otz;Z) 


(55.1) 

(55.2) 


P  =  l  +  eu  (- £  +  1^4-  (("fe)*' +  M9‘)  1“  + 

-u  4-  0(«.‘) 

/>  =  Vfe.'.  +  - ( 1  - Mo.'f.  - c "feVdWo; )3.*)}5i 

- ?  *((  ►feT+^Oy.  +**)/&  +  Va  - 

(£(8j2"'  X<30  f*‘  +  °C«‘o)  (55.4) 

These  variables  in  the  limit  as  tend  to  infinity  must 

match  the  variables  obtained  from  the  equations  away  from 
the  axis  as  the  variable s^jJ^  vanish.  If  matching  is  ac¬ 
hieved,  then  we  have  proof  that  the  expansion  for  the  region 
near  the  axis  was  the  correct  one.  Since  in  this  analysis 
a  composite  expansion  for  the  near  and  far  away  from  the  axis 
regions  is  not  sought,  then  a  matching  of  the  solutions  for 
a  special  case  will  be  adequate  to  indicate  the  correctness 
of  the  expansions  and  will  be  easier  to  perform  than  a  gen¬ 
eral  matching  procedure.  The  special  case  chosen  is  that 
for  a  spherical  shock.  For  this  case  a*  -  O  and  'the  inner 
expansions  can  be  written  as 


f  -  v  iu ^ 


(56.1) 


v  =  6i(i  -  4  0(c? 


'*  (56.2) 
'  (56.3) 
(56.4) 
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For  the  away  from  axis  region,  we  have  the  following  equa¬ 
tions  for  small  9,  9* 

P  =  1  '  +  0fe“)  (57.1) 

f  =(l  -  +  =>C€j2)  (57.2) 

U  =  6“'A  OfeJ  (57.3) 

V  =  CtoCSk/^)1  +•  0(€i'ey*)  (57.4) 


Following  the  same  matching  procedure  as  that  for  the  exter¬ 
nal  layer  flow,  we  define  a  variable^  which  describes  an  in¬ 
termediate  region  where  both  the  away  from  axis  and  near  axis 
expansions  are  valid. 

Since  9  r  oCl)  and  9  *  OfeiJ*)  then  <  0(1) 

In  this  region,  we  have 

9  =  f  «  nqn/ecl1 

then  llkv)  9  -♦  O  and  Jiwi  C4  -*  oo 

Applying  the  above  limiting  procedure  in  the  region  of  common 
validity,  we  must  have 

^  C  -  ft )  o 

€t'o  ■+  O 

>/here  £o  and  § C  denotes  the  various  expansions  for  the  vari¬ 
ables  for  the  outer  and  inner  region  respectively. 


66 


Matching  v 


In  the- intermediate  region,  it  is  known  that 

F  =  ~  Mf) 

Therefore  the  above  is  zero  in  the  limit,  assuring  matching. 


Matching  U  we  can  write 


Again,  making  use  of  the  expression  for  <p>  we  obtain  match¬ 
ing  for  U.  . 

Matching  p 


Introducing  again,  we  then  achieve  matching  for  p 

The  above  matching  process  has  demonstrated  the  correctness 
of  the  near  axis  expansion.  ':Ie  now  use  the  more  general  ex¬ 
pansions,  (55.1  to  55. A)  to  determine  the  layer  thickness  and 
contact  surface  pressure  distribution  near  the  axis.  The 
results  are  ■ 
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tj.  -  € ( v/l6q;HiKqgt> 

+  &/z)z/B*i  “  i^/a-(i(c-(c'- 1)^)  -  (0LfiV2 


-  £;6(l-aaO0/Qi  ((ca-i)'*-  c.) 

+g)a(c-u«-iy*)  + 


C^iViX^-Cc*-!)’^)  <■  (- t^jaQ^,, (-|-(cc'. if 
-C3)  -  Cc*-lV/j)  -B/,)  -  6,3^ 

(  ***  ”  aai)  ~  *  fk)\ $&$(¥■ c " ' (c ' "  ^  ~ 

Cc*-  v)^))  -^i*+J-)S(||^l(ic»-iCc*-i)5A)  -£(c*  -Cc-^H))^*- 


P"  1  *  *<•(*<*■  +  VU4«a-  $^C<P^+M«0#.  +  2Qi)‘ 

#6  *  -  H»t  *  +  H("Sk,T  ■*•  *^>0  B  =/‘t(''te)*  +  4a.-) 

+  4<u) -  '  ^  ' 

e*  =  *  hfe-y  +  Mg;  t  = 

"{g=f +  ^ai 


These  expressions  (  ^  and  p  )  will  then  be  used  in 
Appendix  (A)  to  derive  matching  conditions  for  the  internal 
and  external  layer  near  the  axis  of  symmetry.  From  this  pro¬ 
cedure,  we  can  extract  internal  shock  and  external  (shock  ra- 
dii  of  curvature  and  the  positions  of  the  shocks  in  space. 

This  information  //ill  be  sufficient  starting  conditions  for 
an  integration  across  the  layer  of  the  away  from  the  axis 
equations.  By  this  procedure,  the  full  bow  region,  from  the 
axis  to  the  corner  region,  may  be  calculated.  • 
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Discussion  OF  DDSULTS 


From  the  analytical  bow  region  results  scaling  parameters  are 
defined  which  can  be  used  in  experimental  simulations  and  in 
extensions  of  calculated  results  to  other  conditions.  They 
can  be 'obtained  by  examining  the  first  order  results  of  the 
bow  analysis.  For  the  internal  flow,  we  find  to  the  first  ap¬ 
proximation  that  the  geometry  of  the  layer  (i.e.,  shock  posi¬ 
tion  and  layer  thickness)  depends  on  the  parameters  h  ,  £;0 
j)  and  Rsio  (TT,©co  ,  n  ,  .  To  this  list, 

we  must  add  external  layer  parameters  and  RSe  .  How¬ 

ever,  from  the  matching  of  contact  surface  position,  it  is 
determined  that  R9e  is  a  function  of  RSt'#J©eo,n)€i0  and  there¬ 
fore  only  €e  is  added,  giving  as  the  significant  defining  pa¬ 
rameters  of  the  bow  region  ©00,  r>,  ,  €e  t  RSl;0  •  These 

may  be  further  reduced  if  we  consider  very  high  Ilach  number 
flow  for  which  and  €«»  reduce  to  &j*  and  .  The  geo¬ 
metry  of  the  flow  is  now  independent  of  stream  liach  number. 

For  a  correct  experimental  simulation,  all  of  the  above  para¬ 
meters  should  be  matched  between  the  actual  system  and  the  ex¬ 
perimental  setup. 

Also,  for  a  given  system  with  fixed  ambient  and  exhaust  gas 
composition  (i.e.,  fixed  )  and  nozzle  conditions  (i.e., 

fixed  ©»,n  )  the  bow  geometry  will  scale  with  RSt;0  or  Tf1/afor 
all  flight  conditions  (i.e.,  altitude,  velocity)  and  system 
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thrusts  (i.e.,Foj).  This  last  effect  is  illustrated  in  figure 
4  which  compares  in  a  coordinate  system  which  is  reduced  by 
the  scale  Rsi0  experimentally  determined  bow  region 

geometries  obtained  at  different  values  of  rr .  Also  shown 
in  figure  (4)  is  a  comparison  between  experimental  results 
carried  out  at  nearly  the  same  tr  but  with  different  exter¬ 
nal  flow  Mach  numbers.  As  can  be  seen,  the  agreement  in 
positions  both  for  varying  TT  and  M«e  is  very  good,  indi¬ 
cating  that  the  scaling  parameters  are  correct,  and  also  that 
for  high  Mach  numbers  the  geometry  is  effectively  independ¬ 
ent  of  stream  Mach  numbers. 

The  correctness  of  the  scaling  is  further  substantiated  by  ex- 
amining  figure  (6)  in  Charwat  and  Faulmann  .  In  this  figure, 
many  experimental  values  of  axial  bow  layer  thickness  divided 
by  ^sio  are  plotted  versus  Tt  .  For  both  =  2.75  and  7.1 

it  is  shown  that  the  above  ratio  becomes  constant  as  Tr  becomes 
large,  indicating  the  correctness  of  the  scaling  parameter 
Rsio  .  It  should  also  be  mentioned  that  the  value  of  the 
ratio  reached  forTY»l  and  M<ec  =  7.1  is  in  agreement  v/ith 
the  value  obtained  from  figure  (4). 

Having  established  that  the  geometries  of  the  shocks  and  con¬ 
tact  surface  are  fixed  in  a  reduced  coordinate  system,  we  then 
note  that  to  first  approximation  all  the  nondimensionalized 
flow  variables  will  also  scale  since  they  are  only  functions 
of  <p.>  R/ Rjto  )  ^ >  © oo;  i  wh  i  c  h  are  fixed  quantities  for  any 
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given  point  in  the  flow  if  the  nozzle  conditions  and  exhaust 
and  ambient  gas  composition  are  fixed. 

Another  result  of  the  analytic  analysis  of  the  bow  region  is 
that  it  defines  the  accuracy  of  the  Newtonian  impact  approxi¬ 
mation  which  was  applied  by  Laurmann6  to  this  problem.  In 
his  technique,  the  contact  surface  position  is  determined  by 
the  balancing  of  Newtonian  impact  pressure  along  the  contact 
surface . 

This  approximation  has  been  shown  to  yield  good  results  for 
surface  pressure  for  the  external  bow  region  flow  25 . 

For  the  internal  bow  region  layer  flow,  we  must  examine  our 
analytic  solution  to  establish  the  accuracy  of  the  Newtonian 
impact  theory  in  describing  the  surface  pressure.  From  the 
first  order  results,  the  layer  thickness  is  a  constant.  Also, 
to  first  order  p-essure  is  constant  across  the  layer.  The 
actual  contact  surface  pressure  will  then  be  determined  by  the 
shock  shape,  which  for  constant  layer  thickness  is  'the  same 
as  the  contact  surface.  Therefore,  we  vrould  expect  the  New¬ 
tonian  impact  analysis  to  give  reasonably  good  results  for 
the  internal  as  well  as  the  external  flow  contact  surface  pres' 
sure.  The  contact  surface  shape  in  the  bow  region, as  deter¬ 
mined  by  the  matching  of  the  impact  pressure  aclross  it,  should 
then  be  reasonably  well  predicted  by  Laurmann’s6  analysis. 
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To  test  this  hypothesis  the  contact  surface  position  predict¬ 
ed  by  Newtonian  impact  analysis  is  compared  in  figure  (5)  with 
the  predictions  made  in  this  study  and  in  reference  (7)  for  a 
characteristic  oow  interaction.  The  agreement  shown  is  very 
good,  indicating  the  usefulness  of  this  approximation  to  quick¬ 
ly  obtain  bow  region  geometry.  The  success  of  the  technique 
in  this  axisymmetric  case  gives  one  confidence  in  extending 
this  simple  technique  to  asymmetric  bow  region  flows  where 
the  angle  of  attack  of  the  jet  to  the  free  stream  is  not  large. 
For  the  transonic  corner  and  supersonic  far  field  regions, 
Laurmann’s  6  technique  must  necessarily  fail  due  to  the  tnick- 
ening  of  the  layers  and  the  acceleration  of  the  flow  in  these 
regions. 

A  number  of  calculations  of  the  bow  region  flow  have  been  car¬ 
ried  out,  utilizing  the  analyses  outlined  in  the  preceding 
section  and  in  Appendix  A.  These  calculations  have  attempted 
to:  (i)  compare  predicted  bow  interactions  with  results  from 
an  existing  numerical  technique;  ( ii )  compare  predictions 
with  an  experimental  result  by  Zakkay7;  (iii)  determine  the 
variation  of  bow  region  geometries  and  pressure  distributions 
with  variations  in  ©<*  and  *j  (primary  jet  flow  parameters); 
(iv)  give  detailed  predicted  properties  for  two  cases  which 
are  characteristic  of  "cold"  jet  and  actual  jet  operation. 


In  the  comparison  of  the  present  predictions  with  those  of 
Rudman  and  Vaglio-Laurin^  identical  stream  and  undisturbed 
plume  parameters  were  utilized.  Predictions,  which  included 
terms  with  errors  0(€te)  and  o(£t'©Jwere  made  of  flow  field 
geometry  and  pressure,  as  well  as  the  other  variables  of  in¬ 
terest.  As  can  be  seen  from  figure  (6),  the  present  results 

1- 

and  those  calculated  by  the  numerical  technique  of  reference 
(7)  do  not  agree  very  well  in  position,  although  shapes  are 
similar.  Also,  the  comparison  of  pressure  distributions  in 
figure  (7)  do  not  agree  very  well.  There  are  two  possible 
sources  of  error  which  could  cause  the  poor  agreement  shown. 
The  first  possibility  is  that  the  example  calculated,  with 

=  1.4  is  just  too  difficult  a  test  for  the  present  the¬ 
ory,  which  is  based  on  having  Be,#,- ft  1.0.  Since  the  internal 
flow  solution  proceeds  in  half  powers  of  Gie  then  for  the  ex- 
ample  calculated  terms  higher  than  the  first  are  really  not 
very  small.  This  causes  poor  convergence  of  the  solution, 
which  is  indicated  by  diminishing  oscillations  of  the  vari¬ 
ables  about  their  correct  values  as  higher  order  terms  are 
added  to  the  solution.  This  type  of  behavior  is  observed  in 
figure  (7).  However,  figure  (6)  does  not  show  this  effect; 
and,  therefore,  this  cannot  be  the  total  explanation  for  the 
lack  of  agreement. 

A  possible  source  of  error  in  the  technique  of  reference  (7) 
is  the  imposition  of  a  point  interaction  (i.e. ,  non  penetra— 


tion)  for  the  corner  region  in  the  calculations  of  reference 
(7).  As  will  be  explained  further  in  the  section  on  the  cor¬ 
ner  region,  the  point  interaction  assumption  puts  a  constraint 
on  the  bow  region  calculation,  which  could  possibly  over— em¬ 
phasize  the  corner  region  flows  upstream  influence  on  the  bow 
layer  properties.  This  could  result  in  the  large  difference  in 
internal  layer  thickness  which  is  observed  in  figure  (6).  Tak¬ 
ing  the  above  two  effects  together,  it  is  possible,  within  the 
bounds  of  error  of  the  present  results,  to  show  agreement  bet¬ 
ween  geometry  and  pressure  distribution. 

It  should  be  mentioned  that  the  results  of  reference  (7)  have 
been  compared  favorably  with  the  experimental  results  of  Zak- 
kay  .  However,  the  comparison  is  not  definitive  since  the 
analytical  plume  flow  model  used  had  an  exponent  of  2.5,  which 
has  been  shown  in  the  section  on  jet  exhaust  models  to  be  in 
large  error  for  the  case  calculated. 


Good  agreement  between  the  prediction  of  bow  region  geometry 
and  the  experimental  results  of  Zakkay7  is  shown  in  figure 
(3).  No  linear  scale  was  given  in  the  experimental  schlieren 
photograph;  and,  therefore,  it  was  necessary  to  arbitrarily 
scale  the  experimental  and  predicted  results  by  setting  the 
jet  to  internal  shock  distances  equal  to  each  other.  The 
contact  surface  pressure  distribution  for  this  same  case  is 
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shown  in  figure  (9).  Two  other  calculations  we re  made  for 
these  same  conditions  with  the  only  difference  being  that 
in  the  analytic  plume  model}  exponents  H  =  5.0  and  2.5  were 
used  instead  of  the  more  correct  4.17  value.  In  figure  (10), 
these  -predictions  show  a  marked  sensitivity  of  bow  region  geo¬ 
metry  (i.e.,  shock  and  contact  surface  position)  to  plume 
model  exponent.  The  variation  of  the  bow  geometry  with  plume 
exponent  is  larger  than  the  error  bounds  of  the  solution  in¬ 
dicating  that  if  accurate  predictions  of  bow  region  properties 
are  to  be  obtained  then  the  plume  model  used  must  also  be  very 
accurate . 

In  noting  the  good  agreement  between  the  predicted  and  experi¬ 
mental  geometry,  one  is  tempted  to  state  that  the  calculated 
flow  properties  must  also  be  in  good  agreement  with  experi¬ 
mental  values.  However,  this  is  risky  in  that  experience  with 
co-flowing  plumes  and  the  good  agreement  of  the  calculated 
results  of  reference  (7),  using  an  incorrect  plume  model,  with 
experimental  results  show  that  the  coarse  geometry  of  the 
interaction  (i.e.,  shocks  and  contact  surface)  can  be  pre¬ 
dicted  relatively  easily  but  in  all  probability  the  detailed 
experimental  flow  properties  are  not  so  easily  matched  by  cal¬ 
culations  . 

Figures  (11),  (12),  (13)  and  (14)  show  the  effect  of  varying 
©00  and  fc’j  on  bow  geometry  and  contact  surface  pressure 
distribution. 
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From  figure  (11)  it  can  be  seen  that  increasing  ©eo  f rom  125.5 
degrees  to  185.0  degrees,  which  could  be  accomplished  by  eith¬ 
er  increasing  ©n  or  decreasing  Mj  ,  results  in  a  definite 
increase  in  layer  thickness  both  for  the  internal  and  external 
flow.  ■  Since  external  layer  thickness  is  directly  proportion¬ 
al  to  contact  surface  radius  of  curvature,  then  an  increase  of 
layer  thickness  indi :ates  an  increase  in  radius  of  curvature. 
Consistent  with  the  increase  in  contact  surface  radius  of  cur¬ 
vature,  we  note  from  figure  (12)  an  increase  in  contact  sur¬ 
face  pressure  with  an  increase  in  ©«.  In  general  increases 
in  ©«  with  all  other  parameters  fixed  give  increases  in 
contact  surface  radius  of  curvature,  yielding  a  more  blunt  bow 
region  interaction.  In  figures  (13)  and  (14)  we  can  see  the 
effect  of  increasing  Yj  with  all  other  parameters  fixed.  In¬ 
creasing  £j  from  1.2  to  1.4  in  figure  (13)  causes  a  definite 
thickening  of  the  internal  and  external  layers.  Also,  from 
the  contact  surface  pressure  curves  in  figure  (14)  we  note 
that  increasing  results  in  an  increased  pressure  level. 

Both  of  these  effects  are  primarily  the  result  of  an  increase 
in  radius  of  curvature  of  the  contact  surface.  In  conclusion, 
an  increase  in  either  ©coor  for  all  other  parameters  fixed 
causes  an  increase  in  radius  of  curvature  of  the  contact  sur¬ 
face  resulting  in  a  more  blunt  bow  interaction.  In  this  case, 
pressure,  density  and  temperature  within  the  layers  do  not 
decay  as  rapidly  as  in  the  less  blunt  interaction  case. 
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In  figures  (15-13)  detailed  distributions  of  flow  properties 

for  a  "cold*'  jet  and  actual  jet  operation  are  given.  The  cold 

/ 

jet  case,  shown  in  figures  ( 15,16)  would  be  representative  of 
a  wind  tunnel  simulation  where  the  ratio  of  jet  to  freestream 
stagnation  temperature  is  2.0  and  air  is  used  as  the  jet  gas. 
The  actual  jet  case  shown  in  figures  (17,13)  is  representative 
of  a  real  system  that  is  operating  at  150,000  ft.  altitude, 
whose  jet  gases  are  composed  primarily  of  H20  and  C02  which 
is  modeled  by  setting  =  1.24  .  The  ratio  of  jet  to  free¬ 
stream  stagnation  temperature  is  .75,  which  is  representative 
of  the  actual  system. 

As  can  be  seen  from  figures  (15)  and  (17),  the  flow  geometry, 
streamlines  and  isobars  for  the  two  cases  are  very  similar. 
However,  from  figures  (16)  and  (13)  we  note  wide  differences 
between  the  two  cases.  In  figure  (16)  it  can  be  seen  that 
the  velocity  discontinuity  between  the  internal  and  external 
flow  is  very  small  at  the  contact  surface  and  consequently, 
mixing  effects  along  it  will  be  small.  \Je  also  note  a  mode¬ 
rate  jump  in  density  as  we  go  from  the  external  to  internal 
flow.  Since  pressure  is  equal  across  this  surface,  then  the 
above  result  indicates  that  we  have  a  less  dense  and  hotter 
external  layer  flowing  past  a  more  dense  and  cooler  internal 
layer.  In  figure  (13)  we  note  some  significant  differences 
from  the  "cold"  jet  case  of  figure  (16).  Here,  we  find  that 
there  is  a  large  velocity  discontinuity  which  leads  to  strong 


77 


nvnijifipB! 


'Wf!" 5f -t4 hi IM1* 4-J-L u kji( ■ 1 1? .^wwi rV’*^-5 


mixing  effects.  Also,  the  very  large  discontinuity  in  den¬ 
sity  indicates  that  a  very  hot  and  low  density  external  layer 
is  flowing  past  a  much  cooler  and  denser  internal  layer. 

If  heat  transfer  were  allowed,  there  would  be  a  heating  up  of 
the  cooler  internal  layer  by  the  hot  external  layer.  It  can 
be  concluded  that  in  the  actual  case  mixing  and  heat  trans¬ 
fer  effects  will  be  much  larger  than  the  experimental  simula¬ 
tion  case. 
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ANALYSIS  0?  CQ3NE3  REGION 

For  an  exact  treatment  of  this  region,  the  inviscid  external 
and  internal  layer  bow  flows,  as  well  as  the  relativel17'  mo¬ 
tionless  corp  and  the  viscous  mixing  region  separating  the 
core  and  bow  flows,  must  be  simultaneously  calculated.  nhe 
reasons  for  utilizing  this  coupled  approach  have  been  discus¬ 
sed  previously  in  references  (2)  and  (3).  Such  a  complex 
treatment  is  beyond  the  scope  of  the  present  study  and  sim— 
plif ied  analyses  are  applied.  Unlike  the  bow  region  and  down¬ 
stream  region  analysis  (to  be  discussed),  the  corner  region 
flow  is  not  amenable  to  analytic  treatment  because  of  its 
transonic  nature  and  complicated  coupling  with  the  core  flow 
and  the  viscous  layer  separating  these  regions. 

An  approach  is  then  taken  which  will  provide  an  estimate  of 
oome  of  the  flow  characteristics  within  the  corner  region  and 
will  also  assess  the  significance  of  the  corner  region  flow 
on  the  bow  and  downstream  region  solutions.  In  brief,  this 
section  will  outline j  (i)  the  significance  of  the  corner 
region  flow  on  the  bow  region  and  downstream  flow  properties 
as  calculated  by  their  respective  analyses;  (ii)  an  engineer¬ 
ing  estimate  of  the  pressure  level  in  the  core  region;  (iii) 
the  location  of  the  plume  boundary  based  on  this  core  pres¬ 
sure;  (iv)  the  qualitative  penetration  of  the  intercepting 
shocx  layer  mass  flow  into  the  internal  layer,  and  the  loca¬ 
tion  of  the  internal  layer  sonic  line. 


Before  establishing  the  significance  of  the  corner  region  flow 
to  the  bow  and  downstream  flow  solutions,  a  brief  description 
of  the  flow  in  the  corner  region  will  be  given.  Referring  to 
figure  (19),  the  jet  gases  upon  exiting  the  nozzle  expand  in¬ 
to  the'  low,  nearly  constant  pressure  core  flow.  An  intercept¬ 
ing  shock  and  shock  layer  form  which  are  identical  to  those 
for  a  plume  in  a  quiescent  ambient  of  the  same  pressure.  This 
intercepting  shock  and  layer  eventually  intersect  the  strong 
internal  shock,  resulting  in  a  triple  point  shock  configuration 
similar  to  that  found  at  the  Mach  disc  for  plumes  in  a  quies¬ 
cent  or  co-flowing  ambient.  A  complex  system  of  shocks  and 
expansions  are  then  generated  within  the  internal  and  inter¬ 
cepting  shock  layer,  which  accelerates  and  turns  the  flow  in 
the  downstream  direction.  The  viscous  layer  separating  the 
core  and  internal  flow  reattaches  at  a  downstream  location, 
recompressing  the  flow  to  a  pressure  greater  than  that  of  the 
core.  All  of  the  above  flow  processes  interact,  establish¬ 
ing  a  unique  geometry  for  the  corner  region.  ' 

Having  defined  the  corner  interaction  an  order  of  magnitude 
assessment  of  its  effect  on  the  bow  and  downstream  flow  so¬ 
lutions  will  now  be  given. 
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Significance  of  the  Corner  Region  cn  the  Bov;  and  Downstream 
Solution  “  - 

Since  the  bow  region  layers  are  thin  and  bounded  by  strong 
shocks,  then  the  flow  properties  within  the  layers  are  stronr— 
ly  influenced  by  the  local  shock  conditions  and  only  weakly 
influenced  by  the  conditions  in  the  corner  region  near  the 
sonic  point.  V/e  can  therefore  conclude  that  the  bow  layer 
properties  can  be  calculated  independently  of  the  corner  re¬ 
gion  without  incurring  serious  error.  Since,  as  is  discus¬ 
sed  in  the  next  section,  the  downstream  flov;  depends  on  the 
balancing  of  jet  and  external  flow  interaction  forces  along  the 
contact  surface,  then  the  significance  of  the  corner  region 
flow  on  the  downstream  solution  is  assessed  by  determining  the 
axial  force  balance  that  would  occur  if  the  corner  region  flow 
is  included  or  neglected.  For  the  case  when  it  is  neglected, 
the  jet  flow  is  assumed  to  expand  to  its  vacuum  limiting  angle, 
Oco  . 

The  particular  example  studied  herein  is  that  given  in  a  report 
by  Rudman  and  Vaglic-Laurin^.  The  axial  force  created  by  the 
jet  flov;  is  given  by 

F*  -  -  jf  %  ^  V  ^  ^  -  jpn*cAA 

Applying  this  to  the  cases  with  and  without  the  corner  region 
and  taking  their  difference  divided  by  the  total  force  of  the 
jet,  v;e  obtain  e« 

aFx  =.  ~  2  '"9 

f  «s"ms wycwyJ? 

2.Gco 
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where  *  ±  J  <  1  and 

Substituting  in  the  numerical  values  Ooo  = 
grees,  n  =  4.17,^c°3^  =  ^  we  obtain 


0  <  »  <  9rs 

130. 45>  =  55  de- 


=  .0425  or  4  percent 


AFx 

Fx  t«t*.i 

This  error  is  well  below  that  created  by  using  the  approximate 
analytical  source  flow  expression  (l)  which  in  some  cases  can 
be  in  considerable  error 

Since  the  details  of  the  corner  region  flow  are  not  signifi¬ 
cant  to  the  calculation  of  bow  or  downstream  properties,  then 
engineering  estimates  of  core  pressure,  plume  boundary  and 
penetration  will  be  sufficient  to  approximate  the  corner  re¬ 
gion  flow  for  this  study. 

Core  Pressure 

Qualitatively,  the  external  and  internal  flows  do  not  bend  im¬ 
mediately  downstream  after  passing  the  intercepting  shock  la¬ 
yer  boundary  due  to  the  considerable  lateral  momentum  impart¬ 
ed  to  the  layers  in  the  bow  region.  Instead,  the  external 
shock  moves  continuously  outward  along  with  the  bulk  of  the 
external  and  internal  layer  mass  flow,  leaving  behind  in  its 
"wake"  a  low  pressure  (with  respect  to  the  shock  value)  flow 
adjacent  to  the  core  region.  On  this  basis,  we  would  expect 

the  core  pressure  to  be  low  and  of  the  same  order  as  ambient 
pressure. 
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In  determining  the  core  pressure,  the  experimental  results  of 
Finley  ;  Charwat  and  Faulmann  ,  and  Jarvinen  and  Adams  will 
be  utilized.  From  figure  (7)  of  Finley,  it  can  be  seen  that 
as  ^  as  increased,  the  core  pressure  drops  to  a  nearly  con¬ 
stant  value.  This  behavior  is  also  consistent  with  figure  (11) 
in  Jarvinen  and  Adams  and  figure  (15)  in  Charwat  and  Faulmann. 
In  our  analysis,  which  is  the  limiting  case  when  It  L ,  the 
core  pressure  will  then  be  at  a  constant  level  for  all  values 

of  rr  . 

These  results  are  for  exit  Ilach  numbers  of  1.0,  2.6,  3.1,  3.9, 
and  4.3  and  freestream  Ilach  numbers  of  2.5,  2.75,  .6  and  2.0 
V/e  therefore  conclude  that  this  behavior  with  TT  does  not  de¬ 
pend  on  exit  or  freestream  Mach  number. 

Noting  the  lack  of  dependency  of  the  core  pressure  on  Tf  for 
Tt»l  we  can  now  correlate  all  experimental  results  to  estab¬ 
lish  the  best  value  for  P  core.  From  figure  (20),  it  is  seen 
that  the  most  general  result,  especially  for  high  external 
Mach  number,  would  be  Poore  =  1.  This  result  is  reason¬ 
able  considering  the  wide  variety  of  under  which 

p 

— -&QJ2S.  has  been  obtained  and  also  in  consideration  of  the 
nonuniformity  in  pressure  throughout  the  core  region. 


Plume  Boundary 

As  v;as  mentioned  previously,  the  plume  develops  in  a  re,  -- 
of  nearly  constant  pressure  and,  therefore,  it  is  analogous 
to  the  case  of  an  underexpanded  jet  into  a  still  ambient.  Va¬ 
rious  analytical  approaches  describing  the  location  of  the 
boundary  of  the  plume  have  been  developed  .  (for  eg.  see  re¬ 
ference  2?  and  reference  2f),  Host  of  them  give  reasonable 
results  near  the  nozzle  exit  but  diverge  considerably  from 
the  actual  boundary  away  from  the  exit.  For  this  study,  an 
analytical  technique  '  is  used  which  assumes  a  fixed  plume 
shape  for  all  nozzles  when  the  coordinates  are  nondimension- 
alized  with  respect  to  plume  boundary  maximum  radius  and  its 
axial  coordinate.  The  advantage  of  this  method  is  that  it 
gives  the  correct  boundary  location  not  just  near  the  nozzle 
exit  but  far  from  it  as  well.  The  method  relies  on  setting 
up  momentum  and  force  balances  between  the  exhaust  and  core 
region  gases.  The  shape  of  the  plume  when  nondimensionalized 
is  represented  well  by 

ty/#  ”  ■  ? J  ■  ((X/ ~  V^n  (X/Xmy  J 
* 

v/here  j  =2.5  gives  the  "best"  fit  with  numerical  data.  The 
balance  of  lateral  and  axial  momentum  give  respectively 

Peoce J  ~  Poj  A# 

Pcore  TTT^Xyn  ~  Foj  A*  f  2 


where  fi  end  Jx  are  functions  of  the  nozzle  geometry  and  ex— 
haust  gas  conditions.  From  the  above  expressions  and 

can  be  found  for  any  jet  exit  and  ambient  condition, 
thereby  giving  the  plume  boundary.  This  technique  was  uti¬ 
lized  in  locating  the  plume  boundary  for  the  example  given 
m  Ptudman  and  Vaglio— Laurin' s  study.  Cood  agreement  is 
found  with  experimental  results  as  can  be  seen  in  figure  (3). 
Since  we  are  dealing  with  highly  underexpanded  plumes  (i.e., 

)  where  the  intercepting  shock  layer  is  thin,  then  this 
boundary  location  is  also  assumed  to  be  a  good  approximation 
for  the  intercepting  shock  location. 

Penetration  of  the  Internal  Laver 

Referring  to  figure  (19),  We  can  establish  a  qualitative  des¬ 
cription  of  the  corner  region  flow.  The  reflected  shock  ema- 
nating  from  the  triple  point,  at  a,  crosses  the  intercepting 
shock  layer  and  intersects  the  plume  boundary  at  c.  Since 
the  boundary  must  be  at  constant  core  pressure,  then  the  shock 
must  be  reversed  and  reflected  as  an  expansion  c-b.  Crossing 
the  layer,  this  expansion  then  intersects  the  slip  stream  sur¬ 
face  at  b,  being  once  again  reflected  and  partially  trans¬ 
mitted  as  an  expansion  into  the  slightly  supersonic  internal 
layer  flow  d.  The  surface  b-a  remains  nearly  straight  since 
only  a  small  drop  to  sonic  pressure  over  this  region  is  ex¬ 
pected,  which  requires  very  little  turning  along  b-a  of  the 
supersonic  flow  in  region  a-b-c.  The  internal  flow  sonic 
point  occurs  at  b  because  if  it  were  sonic  before  this  point 
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the  non  divergence  of  the  streamlines  in  this  region  would 
unrealistically  prevent  the  flow  from  accelerating  further. 
Also,  the  sonic  point  cannot  lie  beyond  b  because  the  di— 
verging  nature  of  the  flow  would  not  allow  sonic  velocity 
to  be  achieved  from  thfe  initial  subsonic  state. 

./ithin  the  slightly  supersonic  flow  region,  d,  the  trans¬ 
mitted  expansions  b-c  are  reflected  off  the  nearly  constant 
pressure  sonic  line  surface  as  compressions.  These  return 
to  the  intercepting  shock  layer  and  are  transmitted  to  the 
core  region.  The  continual  process  of  reflection,  transmis¬ 
sion  and  interaction  of  these  waves  results  in  the  rapid  turn¬ 
ing  of  the  flow  in  the  downstream  direction  as  is  illustrated 
in  figure  (19).  Some  idea  of  the  shape  of  the  sonic  line  and 
its  limiting  characteristics  can  be  obtained  by  applying  the 
ideas  developed  by  Hayes  and  Probstein  for  blunt  body  flows. 

Taking  od  as  the  angle  that  the  contact  or  slip  stream  sur¬ 
face  makes  with  the  sonic  line  at  these  respective  surfaces, 
we  can  write  ' 


is  streamline  velocity,  vorticity  and  S  and  n 

are  the  coordinates  along  and  perpendicular  to  the  slip  stream 
and  contact  surface.  At  point  b  ve  have  nearly  a  straight 
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slip  stream  b-a  and  therefore  the  second  term  in  to  - 
negligible.  Also,  the  shock  is  nearly  spherical  and  in  the 
vicinity  of  point  a  the  vorticity  generated  is  small.  Under 
these  conditions,  (Mi- °  and  therefore  cO  ^  at 

b.  At  the  contact  surface,  point  e  ,  different  conditions 
prevail.  Using  the  first  order  expressions  derived  in  the 
bow  region  analysis  section  and  assuming  for  simplicity  that 
the  internal  shock  and  contact  surface  are  spherical  in  nat¬ 
ure,  then  we  can  write  in  terms  of  nondimensionalized  variables. 

%  *  "(ikf?  %.  - 

7 

Using  Rudman's  example,  Wo  can  write 

also,  we  have 


^  geo=  1.54 


The  ratio  of  these  expressions  gives 

tan  =  1.^5  which  for  <p  =  55°  gives  u)  =  62  degrees 

Now  that  the  sonic  line  angles  are  estimated  at  both  bound¬ 
aries,  we  can  then  sketch  its  location  throughout  the  layc  - 
as  shown  in  figure  (19).  The  initial  turning  of  the  sonic 
line  in  the  downstream  direction  from  point  b,due  to  centri¬ 
fugal  pressure  gradient  effects',  is  typical  of  the  rapid  tran¬ 
sonic  expansion  occurring  in  d  and  gives  the  sonic  line  its 


-u.  i  j«  14  wyu 


1 11. -  -f  *  IT"  <***■ 


characteristic  s  shape.  The  last  limiting  characteristic 
lines  which  originate  in  regions  d  and  f  are  also  shown  in 
figure  (19).  Having  obtained  a  qualitative  estimate  of  the 
penetration  of  the  intercepting  shock  layer  into  the  internal 
layer,  we  can  now  make  an  estimate  of  the  actual  penetration 
for  the  example  calculated  by  Rudman  and  Vaglio-Laurin^.  The 
first  step  is  to  estimate  the  intercepting  shock  layer  thick¬ 
ness  which  can  be  determined  on  the  basis  of  mass  flow  con¬ 
siderations.  Subscripting  layer  properties  by  Jt  and  utiliz¬ 
ing  standard  notation  for  undisturbed  plume  quantities,  we 
can  write 

60s 

pjiUjl  R*  Sm  ©A  Sft  -  j’ j0oot'  Uooi  3.TT R^sih9  dp 

&SL 

Substituting  in  the  analytic  expressions  for  poo:  and  U00L 
which  can  be  derived  from  the  analytical  plume  flow  model, 
and,  letting 


f*i.  ~  4  V2.)  Ub.  =  n°]/z. 


we  obtain  - 

h-  =  a-euLeaa"M»6my 

Rsi°  C05"  §§£  61  f  e&/W7#) 

For  800=130.45°  eA  =  55°  and  6^=1  we 

have  =  .0627 

R  $u 

This  value  is  characteristic  of  plumes  whereTf»l  .  The  se_ 
cond  step  involves  determining  the  reflected  shock  position 
from  the  triple  point  solution.  Using  the  initial  conditions 
Hooi  =  9.2,  >$  =  86  for  the  strong  shock  and  ,5=  20° 
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for  the  intercepting  shock,  we  obtain  the  reflected  shock  po- 
sition  from  D’Attorre  which  is  shown  in  figure  (&). 

Assuming  that  the  Mach  number  downstream  from  the  reflected 

shock  is  constant,  we  can  then  determine  the  initial  reflected 

expansion  wave’s  position  and  intersection  with  the  slip  stream 

surface.  The  position  of  this  wave  as  shown  in  figure  (&) 

defines  the  extent  of  the  penetration  into  the  internal  layer. 

Unlike  the  point  interaction  or  zero  penetration  model  pro- 

7 

posed  by  Rudman  and  Vaglio-Laurin  ,  we  note  a  considerable  pe¬ 
netration  of  the  internal  layer  by  the  intercepting  shock  la¬ 
yer  flow.  This  fact  is  borne  out  experimentally  in  a  study 
on  the  similar  problem  of  a  jet  impinging  normally  and  at  an 

angle  to  a  flat  plate.  The  constraints  on  the  numerical  so- 

7 

lution  of  Rudman  and  Vaglio-Laurin  imposed  by  their  assump¬ 
tion  of  a  point  interaction  (i.e.,  zero  penetration)  for  the 
corner  region  could  be  a  source  of  considerable  error  in  their 
calculation  of  the  bow  region. 

It  should  be  mentioned  that  the  large  penetration  indicated  in 
this  study  does  not  necessarily  mean  that  the  bow  region  flow 
is  strongly  influenced  by  the  corner  region  penetration.  This 
is  because,  in  the  absence  of  the  corner  region  flow,  the  inter¬ 
nal  flow  streamlines  would  penetrate  the  layer  much  like  the 
slip  stream  separating  the  internal  and  intercepting  shock 
layer  flows.  This  can  be  seen  by  comparing  streamlines  ob¬ 
tained  from  the  bow  analysis  in  figure  (17)  with  the  slip 
stream  position  shown  in  figure  (8). 


From  the  results  presented  a  scaling  criteria  for  the  corner 

region  is  now  outlined.  The  major  characteristic  to  be  scaled 

is  the  undisturbed  plume  boundary,  which  for  a  given  nozzle 

and  exhaust  gas  composition  depends  on  P«j/  #  por  an 

exact  simulation  of  the  plume  boundary,  this  ratio  must  be 

matched.  If,  for  a  given  system,  it  can  be  shown  that  the 

plume  boundary  coordinates  are  proportional  to  or  Tf ^ 

then  the  plume  boundary  will  scale  like  the  bow  region.  The 
2Q 

technique  used  to  locate  the  plume  boundary  normalizes  the 

geometry  by  and  the  plume  maximum  radius  and  its 

axial  location.  When  these  are  plotted  as  a  function  of 
j  /  P«»e  it  is  found  that 

Gm.  ,  oc  (Pej/p,e) 

Since  TT  is  proportional  to  Poj/poo.M*^  then 

.Cel  >  x oc  TT  M»e 

rv  r* 

Thus,  if  Moo.  is  fixed,  then  the  plume  boundary  will  scale 
with  Rs,o  .  The  scaling  of  the  bow  region  and  plume  boundary 
then  requires  that  f  ,  ©«»,  H  ,  and  Mcoe  be  fixed  over 
the  entire  flight  range.  Since  does  not  vary  greatly 

over  a  typical  system  trajectory,  then  the  above  scaling  re¬ 
quirements  are  satisfied  for  some  systems  of  interest. 

This  scaling  result  is  substantiated  by  some  experimental’ dl^a 
in  which  the  plume  boundary  triple  point  distance  away  from 
the  axis,  when  divided  by  the  scale  length  RSt'e  is  nearly  con¬ 
stant  for  values  of  IT  from  30  to  627. 
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FAR  FIELD  REGION  ANALYSIS 


The  downstream  flow  region  between  the  external  shock  and  axis 
of  symmetry  consists  of  supersonic  layers  bounded  by  mixing 
regions  and  crossed  by  a  wake  recompression  shock  as  is  illus¬ 
trated  in  schematic  (3) 

.mixing  region 


.mixing  region 

Schematic  (3) 

This  general  problem  cannot  be  handled  by  approximate  analy¬ 
tical  techniques;  therefore,  we  will  seek  to  solve  a  reduced 
problem  within  the  framework  of  the  assumptions  listed  in  the 
introduction.  This  problem  will  involve: 

1)  the  neglect  of  all  mixing  regions  (wake 
and  shear  surface) 

2)  the  neglect  of  the  recompression  shock 


Assumption  1)  will  be  valid  at  moderate  altitudes  and  2)  will 
be  reasonable  since  the  recompre ssion  shock  is  weak  compared 
to  the  external  shock.  A  schematic  of  the  flow  and  notation 


The  numerical  matching  criteria  between  the  internal  and  exter— 
nal  flow  requires  that  static  pressure  and  flow  deflection  must 
be  equal  at  the  contact  surface.  Also,  the  stagnation  pres¬ 
sures  will  be  equal  due  to  the  upstream  matching  condition. 

The  particular  region  of  interest  in  this  downstream  analysis 
will  be  at  a  point  characterized  by  the  distance  L  ,  v;here 
•  i-he  ratio  of  f^9to  L.  will  be  a  significant 
parameter  for  the  downstream  problem  and  we  denote  it  by  the 
symbol  Se  -  Rs/L_  .  Once  again,  as  in  the  bow  region,  we  will 
be  considering  the  case  when  9  6^0j6e« L  An  ap_ 

proximate  analytical  method  of  solution  is  available  for  the 
external  layer  flow  near  the  shock  surface  in  region  1.  This 
method  was  first  applied  by  Cole  33  and  involves  taking  the 

limits  in  the  inviscid  flow  equations  as  the  perturbation  quan¬ 
tities. 

lv]coc  “ CO  ,  £e  O  )  S  e  ->  O 

Mooe^e  "■*  00  Moot  €.e  O(l)  » 

The  solution  found  by  this  procedure  is  valid  to  a  distance  of 
o(CeSeL)  measured  inward  from  the  shock  surface.  For  attached 
shocks,  this  solution  would  be  valid  throughout  the  external 
layer.  However,  if  the  shock  has  a  blunt  nearly  normal  for¬ 
ward  region,  the  external  layer  becomes  thicker  than  the  above 
order  of  magnitude  and  the  solution  breaks  down  away  from  the 
shock  toward  the  contact  surface.  This  is  due  to  the  viola- 


tion  in  the  blunt  region  of  the  assumption  that  the  shock 
makes  small  angles  with  the  free  stream.  As  was  pointed  out 

O  1 

by  Cheng  and  Kirsch  in  the  equivalent  unsteady  case,  a 
region  of  thickness  o(s«0 

,  denoted  2  in  schematic  (A) ,  is 
then  found  near  the  contact  surface  which  requires  a  differ¬ 
ent  expansion  in  order  to  develop  a  valid  solution  in  this  re¬ 
gion.  The  solution  found  in  this  "entropy  wake"  matches  with 
the  outer  layer,  and  the  composite  expansions  for  the  entire 
external  layer  will  be  given. 

For  the  internal  layer,  we  have  two  conditions  imposed  by  the 
external  flow  at  the  contact  surface,  which  can  be  used  to  es¬ 
tablish  the  orders  of  magnitude  of  the  flow  variables  within 
region  3,  the  interior  layer.  They  are: 


Based  on  these  conditions,  we  can  then  prescribe  expansions 
for  the  flow  variables  in  terms  of  perturbation  parameters 
that  will  lead  to  a  system  of  analytical  solutions  for  the 
inner  layer. 


Like  the  external  "entropy  wake"  the  internal  layer  is  of  thick' 
ness  0^5e  .  Utilizing  the  orders  of  magnitude  found  from 

the  contact  surface  and  upstream  bow  region  conditions,  we 


then  develop  expansions  which  will  be  solved  giving  expres¬ 
sions  for  layer  thickness  and  other  variables  of  interest. 
Numerically  Hatching  these  expressions  with  those  of  the  ex¬ 
ternal  flow,  along  the  contact  surface,  we  then  achieve  ana¬ 
lytic  solutions  in  the  far  downstream  region.  These  results 
are  found  to  depend  on  the  upstream  bow  interaction  between 
the  plume  and  external  stream.  Taking  results  from  the  bow 
region  consistent  with  the  accuracy  of  the  downstream  region, 
we  then  find  the  complete  solution  for  the  downstream  flow. 

In  the  following  discussion  for  external  flow,  frequent 
reference  will  be  made  to  the  work  of  Cheng  and  Kirsch 
They  solved  an  unsteady  problem  that  in  many  ways  is  anala- 
gous  to  our  external  layer  problem  by  the  application  of  the 
hypersonic  equivalence  principle.  For  details  of  the  exter¬ 
nal  layer  analysis,  the  Cheng  and  Tirsch  paper  should  be  con¬ 
sulted. 


External  Layer  -  Region  1  Analysis 

Utilizing  the  notation  illustrated  in  schematic  (/*.),  the  fol¬ 
lowing  general  set  of  equations  and  boundary  conditions  for 
the  downstream  problem  in  terms  of  von  Mises  variables 
are  • 


Xe 


'bH'e  Re 

+  \Bs.  - 

Re^Xe 


Ut1  +  Ve‘  +  i*SL-  &L.  =  rk  (  V*) 


tte  -  1  'js  e 


Kf.*  =  f(Ve) 


(1.1)  stream¬ 
line  slope 

(1.2)  continu¬ 
ity 

(1.3)  momentum 
normal  to 
streamline 

(1.4)  Bernoulli 
integral 

(1.5)  entropy 
integral 


Boundary  conditions  at  external  shock 


Pe  *  P®e  p'oe  —  ly  (1.6) 

1  1  Xt  t  -t  '  ' 

fe  =  j^«/ (  fca  +  ^ 1  ’ 

U«  =  U®*(1  “  2.  (Moo*  s»m*vS«  )  (1*3) 

Ve  *  U«e  (  «  <Moo\ -  l) cotfie/fve+O  M«f  ( 1  ’ 9 ) 


Boundary  conditions  at  contact  surface 


Pft  cs 

A 


-  ^ t  cs 

=  d  R.u 


(1.10) 

(1.11) 
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y  - 


•  ’v  •^WT  W^  ’!-,T!p^w  -,vm w  &vf,wj  \Tr$  'i»*™r7r?nw  \nmjrr-fT 


Velocity  parallel  and  normal  to  the  shock  are: 

U  “  L(eo»cos^e 
V  —  CLe  LA  fio  e  s  i  vi  /3  c 


”iw-i .«?;Hrr«Tjra,“  ^ ■  JT;?”  ■  ■v,,i.l>-T’ur. w* w w.ywyy 1 **wr 


(1.12) 

(1.13) 


The  expanded  Rankine-Kugoniot  relations  lead  to  the  follow¬ 
ing  orders  of  magnitude  for  the  flow  parameters  near  the 
shock: 


p« 

O  Cf  «•  Uooe  / 

(2.1) 

/>• 

(2.2) 

Ue 

'<>w  ^  C  lA  eoe  *) 

(2.3) 

/V 

Ve 

'"s“  O  C^e  U.eo») 

(2.4) 

where  €e  is  evaluated  at  the  point  v/here  =  TT/;^  t  These 

conditions  will  be  the  upper  boundary  values  for  the  external 
layer.  From  the  orders  of  the  variables  near  the  shock,  the 
dimension  normal  to  the  shock  can  be  inferred  from  continuity 
considerations . 


"s-  O  (fie  Se.  1-)  (2.5) 

The  dimensions  of  ihe  streamwise  coordinate  are  easily  se^n  to 
be 


x  ~  o(l) 


(2.6) 


- - - - - - -  ■ . . 


■it 


These  orders  of  magnitude,  (2.1  to  2.6),  are  then  used  to  form 
asymptotic  expansions  which  are 

pe  —  PooeUa^  §c  p0e  +■ 

Ue  —  Uaoc  ***  S&*  Uooe  Uoe.  + 

Ve  —  L(ooe  Voe  +  •  • 

tye  =  L  Se  6e  ^»e 

Substituting  these  expansions  into  the  equations  of  motion, 
(1.1  to  1.5),  and  boundary  conditions,  (1.6  to  1.13),  v;e  can 
develop  a  sequence  of  equations  which  can  be  solved  and  then 
summed  to  yield  the  solution  to  any  desired  accuracy  in  terms 

of  Mooe  . 

The  first  order  system  from  the  above  procedure  is 

JL  ,  , 

Rs*  ^Xel  (3.1) 


C  B>e/poe)  =■  O 

oXe 

(3.2) 

po«  R»« 

(3.3) 

v.t  =  ^ 

(3.4) 

The  boundary  conditions  on  the  above  at  the  shock  surface  are: 


£>e  *  * 

"  1 


(4.2) 

(4.3) 


Integrating  (3.1  to  3.4)  and  applying  the  boundary  conditions 
(4.1  to  4.3),  we  obtain: 


^oe 


Voe. 


(5.1) 

(5.2) 

(5.3) 

(5.4) 


Making  the  substitution 


O'  a  i 

=  R../2 

in  the  above,  we  then  obtain  results  identical 


(5.5) 

(5.6) 

to  those  of 


As  previously  mentioned,  these  results  are  expected  to  hold 
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near  the  shock  surface;  however,  their  validity  near  the  con¬ 
tact  surface  has  to  be  examined.  From  the  bow  region  solution 
we  know  that  a  blunt  interaction  occurs.  Consequently,  as 
we  follow  the  shock  surface  upstream  to  the  vicinity  of  the  bov 
region,  we  note  that  the  shock  must  become  nearly  normal  in 
form.  Therefore,  it  invalidates  our  assumption  of  small  free 
stream  to  shock  angle,  initially  made  in  our  above  analysis. 
Under  these  circumstances,  one  would  not  expect  our  expansion 
procedure  to  be  valid  in  such  a  region; and  this  turns  out  to 
be  ti.j  case,  as  was  found  by  Cheng  and  Kirsch  ^  for  the  un¬ 
steady  problem.  What  we  must  do  is  alter  our  expansion  pro¬ 
cedure  so  as  to  be  consistent  with  the  strong  part  of  the  shock 

External  Layer-Region  2  Analysis 

Following  Cheng  and  Xirsch  3/f  we  assume  that  initially/  due 
to  the  strong  bow  interaction  of  the  jet  and  external  flow, 
the  external  shock  moves  out  in  a  blast  wave  manner.  This 
requires  that  the  external  shod:  take  the  form: 

«  Ax«v* 

as  X«-*0  where  A  is  a  function  of  the  upstream  bow  inter¬ 
action.  A  requirement  which  must  be  met  by  the  contact  sur¬ 
face  in  the  vicinity  of  the  blunt  interaction  is  that: 

l’"  Rcs«/f?,.  -*o 
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This  states  that  initially  the  contact  surface  must  move  out¬ 
ward  slov/er  than  the  external  shock,  thus  guaranteeing  the 
dominance  of  the  blast  wave  result  for  the  upstream  interac¬ 
tion.  ■  Up  to  this  point,  we  have  paralleled  the  unsteady 
problem  quite  closely.  The  equivalence  principle  is  most  ac¬ 
curate  near  the  shock  with  the  streamwise  velocity  given  by: 

Ue  lieo«  +“  O  C^®e 

The  error  involved  in  assuming  U*  ~  Uceeis  seen  to  be  small, 
of  error  o(l>loo*§e)  .  However,  as  v/e  approach  the  contact 
surface  we  have  for  the  streamwise  layer  velocity,  from  the 
Bernoulli  Integral. 


U«.  «■ 


The  error  involved  in  taking  U~Ucoe  in  order  to  apply  the 
equivalence  principle  is  seen  to  be  large  for  fce  —  1.  of 

In  order  to  utilize  the  equivalence  principle,  we  then  must 
stipulate  that  &e  cannot  equal  one  but  can  be  near  one  as 
long  as  is  small  enough  so  that  the  term  de  *•  does 
not  become  of  order  one  in  the  downstream  region  of  interest. 
If  these  conditions  are  met,  then  we  can  apply  the  equival- 

3/ 

ence  principle  and,  therefore,  make  use  of  Cheng  and  Kirsch’s 
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results  within  the  entropy  layer. 


Assuming  this  form  01  initial  behavior,  we  have  then  speci¬ 
fied  the  type  of  breakdown  of  the  solution  due  to  the  near¬ 
ly  normal  character  of  the  shock  in  the  vicinity  of  the  nose. 
It  was  f oun^  that  the  expansion  for  ^  based  on  flow  near 
the  shock  for  initial  blast  wave  behavior  would  be  of  the 
form: 

£  ~  o On  Yo*)  4  o  (e,(ln  Yo*)*) 

which  would  then  become  infinite  and  invalid  as  an  expansion 
as  Yo*->0  ,  or  as  we  approach  the  contact  surface.  This 
form  indicates  that  the  expansions  cannot  be  made  uniformly 
valid  by  simply  altering  the  scale  of  Yo*  by  a  power  of  £c  • 
In  his  analysis,  it  was  determined  that  the  scale  normal  to 
the  shock  should  be  in  terms  of  Y°»  to  some  power  which  is  a 
function  of  €e  .  Te  can  then  determine  the  orders  of  magni¬ 
tude  of  the  other  variables  by  applying  entrooy  and  energy 
conservation  principles. 

Ol 

Following  Cheng  and  Kirsch  we  then  establish  a  new  set  of 
variables,  fulfilling  all  the  necessary  conditions  within  the 
region  adjacent  to  the  contact  surface.  These  are: 

Ve  =■  U-coc  Voe  +  *  '  •  (7.^) 
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•  • 


Ue 

P* 

fe 

'&«■ 


LI  floe  +  * 


”  P»e  +■  *  *  * 

~  L  jo  +  •  *  • 


X 

*? 


(7.2) 

(7.3) 

(7.4) 

(7.5) 

(7.6) 

(7.7) 


Substituting  these  variables  into  the  equations  of  motion, 
(1.1  to  1.5),  we  obtain  the  following  set  of  first  order  re¬ 
sults: 


-  o 

& 

(S.i) 

K 

Sj-Oq  C-X,0)  —  C  _ 

(3.2) 

J°°* 

=  pce/a;^ 

(3.3) 

Voe 

ll 

fi 

(3.4) 

The  carrying 

out  of  the  integration  in  (8.1) 

and 

(£.2)  gives 

Fe 

=  ^.Cx) 

(9.1) 

-  ^0e*Cx,O)  a0^/p6e 

(9.2) 

Matching  of  the  inner  to  the  outer  solution  is  carried  out  by- 
determining  the  behavior  of  the  solutions  in  an  intermediate 
region  where  they  are  both  valid.  The  results  ^  of  this 
procedure  are: 
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([%')  *  tI?  ,io  i> 

p...  (#.)*» 

Equations  (10.1)  and  (10.2)  are  nonlinear  second  order  dif¬ 
ferential  equations  for  the  pressure  distribution  and  shod: 
position.  If  v/e  can  obtain  an  inner  layer  solution  which 
gives  RcS  as  a  function  of  pCe  ,  then  we  can  combine  these 
results  into  the  above  to  determine  R.se  and  Res  . 

Internal  Layer  -  Re<~ion  3  Analysis 

Referring  once  again  to  our  flow  diagram  to  establish  the 
system  of  coordinates,  v/e  ’have  the  following  set  of  equations 
and  boundary  conditions  for  the  downstream  internal  layer 

(11.1)  stream 
line  slope 

(11.2)  contin¬ 
uity 

(11.3)  momentum 
normal  to  stream¬ 
line 

(11.4)  Bernoulli 
Integral 

(11.3)  Entropy 
Integral 


problem  in  terms  of  von  Rises  variables 


=  V 

bti 

_ 

Ri 


-  M 


=  o 


Uil  +  Vil  -t-  ^=.  a^c^i) 

<■> 

=  H^c) 
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ternal  shock,  we  find 


pics  s  °  Cfeo£0  UctoL  S*2)  (12.1) 

Since  entropy  is  conserved  along  streamlines,  we  can  es¬ 
tablish  the  order  of  magnitude  of  density  in  the  internal 
layer  near  the  contact  surface  from  (11. 5) 

P  =  P**/p^ 

For  Mooi^l  ,  s\y\x &i  Ct  1  the  shock  conditions,  (11.7) 
and  (11. 6), for  the  streamline  which  wets  the  contact  sur¬ 
face  become 

it*,*) 

V/here  subscript  o  denotes  normal  shock  location  in  bow 
region.  '.Then  substituted  into  the  entropy  equation,  (11. 5), 
the  result  is:  t 

=  O(p*coeo  ^e^/et'o)  (12.2) 

From  the  second  contact  surface  condition,  (11.11),  we  find 

=  d-  Rftcs  =  Ve_  -  o(<5e) 

dx  cLx 
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This  also  requires  that 


Vt/u 


(12.3) 


near  the  contact  surface,  which  is  an  upper  bound  on  this 
ratio.  Another  expression  which  can  be  utilised  is  Bernoulli’s 
equation,  (11.4),  along  a  streamline.  Substituting  (12.1  to 
12.3)  into  this  equation  and  assuming Ui  of  the  form 

Lt  t  —  L^ooi,  +■  U  i 

we  obtain 


IU*  +-  2UoolU(''  +  o(St*U0 °C 

.  _2iU_o(£i°  feoCoUeoi  Sc 

J^o Ota 


Therefore,  we  find  that 


+  2>Se*Uco i  Uc") 
=  Uco^ 


U 


t 


(12.4) 


Therefore,  we  can  conclude  that  for  M.i  to  be  constant  in  the 
first  approximation  in  the  above  flow  system,  it  is  neces¬ 
sary  that  Vi^landJebe  small  enough  so  that  the  term 

Se  0(1)  .  These  are  the  same  conditions  ’which 

are  necessary  for  the  external  layer  solution  to  be  valid 
near  the  contact  surface;  therefore,  we  assume  that  this  will 
be  the  case  in  the  present  problem.  This  leads  to  the  result 

~  o(IAool')  (12.5) 
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From  the  contact  surface  condition,  (12.3),  we  find 

Vi,  -  O  (  Uool  oe)  (12.6) 

V/e  have  no u  established  the  orders  of  magnitude  of  all  inte¬ 
rior  variables  of  interest  near  the  contact  surface.  V.Te  now 
postulate  expansions  for  the  variables  based  on  these  orders 
and  then  substitute  them  into  the  equations  of  motion  to  yield 
a  sequence  of  solutions,  which  can  be  summed  to  give  the  solu¬ 
tion  to  any  order  of  accuracy  for  the  internal  layer  flow. 

The  expansions  are  for  £t*0-*O  0  O  3  M  oio  -V  CO 

Pi  ~  pot  +  *  ‘  * 

Pi.  ~  p«>i’a  * 

vj,  =  U»t  Vot  +■  *  *  * 

Ut  “  Uoot  ■+“**’ 


The  first  order  set  of  equations  is 


^Pot'  -  O 
>  Roc  —  _ 3 _ 

'b'+'T  P°c  R-ot 

feU*vjO  -  o 


(13.1) 

(13.2) 
(13-3) 
(13.4) 
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£<i  Lflb.  L*lL  u.~  kubUtoifexl 


Integrating,  we  then  have 


Pot  -  P°tCx) 

w 

(14.1) 

f?*;  *  s  y  f  d  *4/V 

'  L  J£>  f  oi 

(14.2) 

9°  L  “  pfe  t  i*  /  poj  * 

(14.3) 

Combining  (14.1,  14.2  and  14.3),  we  obtain 


2-  f  Esi*.  i<*>C 


k  .Pot 


(14.4) 


Now  that  we  have  results  for  the  internal  and  external  flow, 
.ve  can  combine  them  to  determine  the  location  of  the  shock 
and  contact  surface,  as  well  as  find  the  values  of  the  pro¬ 
perties  across  the  layers. 


From  the  external  flow,  we  have  the  results: 


P°'«=  (tST  *  (10-2) 

Po.0  (  -  Res1)  =  <*/*.  !10-1) 


Also  from  the  internal  flow,  we  have  the  result 


Since 

dition, 


footUot*.  =  ^ooJoUoot  from  the  Stagnation  point  con- 
then  p0ecs  =  -po^Cx.)  •  "■’•e  can  n°w  combine 
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these  relations  (10.1,  10.2  and  14*4)  from  which  property 
variables  can  be  determined. 

Taking  (10.1)  and  (14.4)  we  have 

£x>o)(rs: 

or 


qb 

2. 


(15.1) 


Combining  (15.1)  v/ith  (10.2),  we  then  obtain  for  the  shock 
shape  a  second  order  non-linear  differential  equation 


For  ^/R.St0  >V  1  ,  which  is  in  the  downstream  region  of 

interest,  (15.2)  admits  a  solution  of  the  form 


C  x 


(15.3) 


When  (15*3)  is  substituted  in  (15.2),  we  find  the  value  of 
C  to  be  given  by 


(15.4) 


Using  (15.4)  in  (10.2)  to  find  po<;Cx)  (14.4)  then  be¬ 
comes 

Vj  ‘*t 

=  ^_  r  p.,v  (15. 

Rc»X  '6  j^t-V 


Combining  'l 5. 4)  and  (15* 5) >  we  tnen  obtain  the  following 
expression  for  the  contact  surface  shape  in  terms  of  the  shock 
shape 


&./( 


cr0 

*4jef 


Jo  ’fbt* 


(15.6) 


V/e  now  have  determined  the  values  of  Rs«  ,  Rcs  in  terms 
of  the  functions  <3^  and  jE£_yjL  .  These 

expressions  are  functions  of  the  blunt  interaction  between 


the  external  and  internal  layers  at  the  extreme  upstream 
position  in  the  bow  region  and  are  a  result  of  the  matching 
of  pressures  across  the  upstream  contact  surface. 


V/e,  therefore,  must  postulate  an  upstream  interact  ion  based 
on  our  solution  for  the  bow  region  flow  which  is  consistent 
with  the  approximations  made  in  this  downstream  region. 

For  the  downstream  analysis,  we  require  €t\,,€e.-*0  Mose-*® • 

These  limits  must  also  be  applied  to  our  bow  region  matching 
conditions  which  are  the  equality  of  pressure  and  flow  deflec¬ 
tion  along  the  contact  surface.  Since  our  bow  region  internal 


lie 


flow  is  a  perturbation  of  the  case  where  the  shock  and  body- 

surface  are  a  spherical  shell,  for  €-c*0  =  0  ,  then  v.re 

make  the  approximation  that  the  upstream  interaction  region 

is  spherical  in  nature.  For  the  external  flow,  it  is  a  well 
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established  fact  that  the  downstream  shock  surface  does  not 


depend  on  the  details  of  the  pressure  distribution  along  the 
contact  surf ace, but  only  on  the  total  integrated  pressure  dis¬ 
tribution  (i.e.,  drag)  in  the  axial  direction,  adding  further 
validity  to  the  use  of  the  above  assumption.  Along  this  sur¬ 
face,  the  pressure  induced  by  the  internal  flow  must  match 


that  of  the  external  flow.  To  first  order,  which  is  the  de¬ 
gree  of  accuracy  for  the  downstream  region,  the  pressure  change 
across  the  shock  layer  is  zero  and,  therefore,  we  can  calcu¬ 
late  the  pressure  on  the  contact  surface  by  simply  finding 
the  pressure  level  behind  a  spherical  shock  surrounding  the  jet 
flow  source.  For  the  case  e  1  ^  Svv^  C  —  1.  » 

the  expression  for  pressure  behind  the  shock  reduces  to 

P<-’  -  -  fcoC  Ucoi  -  pood  PooCe^°Ot 

Poo  to  J 

Non-dimensionalized  with  respect  to  the  external  conditions, 
this  relation  becomes 

.  i 

Pe.fr  t  -  £g>4.  peQ/$ULscu. 

1  peoio  joooc.  U.ooe  *• 

From  the  equality  of  stagnation  pressures,  this  then  reduces  to 


Pcs  i  -  —  pcse. 

paOwo 


From  the  source  conditions,  we  have 


pcs.* 


With  the  pressure  level  and  contact  surface  shape  determined, 
we  have  satisfied  our  matching  criteria  and  our  next  step  is 


to  calculate  the  functions  cra  and 


Mg  • 
J© 


Evaluating  the  integral  from  the  internal  layer  solution,  we 


find 


Vjet  Mjet 

§  d  M^l  -  i  cos^l 

o  J0 

J&  cosv,TI£|_  Jj> 


CD /CjOS 


't'j’et 

dWi  a  j"  <*4-L  = 


In  this  result,  the  upper  integration  limit  requires  some  d"  s- 
cussion.  Strictly  speaking,  in  order  to  calculate  the  pres¬ 
sure  and  density  distributions  over  the  full  range  of  <p 
one  would  have  to  take  account  of  the  effect  of  the  plume 
boundary  ~nd  the  mass  flowing  in  the  plume  boundary  layer  on 
the  internal  shock.  However,  for  highly  underexpanded  jets, 
there  is  relatively  little  mass  flowing  in  this  layer  and, 
therefore,  the  integration  of  the  above  to  the  vacuum  limit, 


—  ©eo 

analysis. 


creates  only  a  small  error  in  the  downstream 


We  now  seek  to  determine  the  value  of  cr0  ,  jt  can  shown 

that  for  <<  l  our  postulated  initial  blast  wave  be- 
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havior  for  the  shock  shape  can  be  written  as 


.'/here  Tfe  to  our  degree  of  approximation  can  be  written  as 


lb  =  ife1  —  3- 


^fre-l  )(te  +  -0 

Tn  the  above  Con  ia  the  nose  drag  coefficient  of  a  hypothe¬ 
tical  blunt  nosed  slender  body,  which  takes  on  the  shape  of 
the  contact  surface.  C  on  is  directly  related  to  the 
blunt  interaction  and  is  written  as 

^  bN  —  ^  d  A  / ^®e  Ucoe 

Since  we  stipulated  that  R  es./^Se  -*  o  as  o 

then  cquafions(lO.l)  and  (10.2)  give 


'A&j  +  4*Rs»  Rsr)Rs»  - 


x/RSl< 


from  which 


<r«  - 


Q  &>N„ 
lb 


>.e  now  calculate  C&n  based  on  our  postulated  upstream  in¬ 
teraction.  For  a  spherical  shell  interaction  with  the  pres¬ 
sure  distribution  given  by 

P  -  J^cce.  Ucbe  Cos 

then 

rT  •  n  x  -  cor.  cp 

-and 


2.@ao 


d  A  *  2.TT  Rsfo1  s^<J  d<p 


•  hich  leads  t 


to 


1  i  i 


•lU^Pi^-g.  »JI!  l !  JJJ  4A  ^  Lpn*-*tu  IP  1  J,W  11.^  .i,.  L .  J,Jj.l|  I  Rv.  J  WP|p|| 


C&n  -  H 


r 


©CD 


cosnn^»Lsi^9v  cos  9*  49* 

2J9cc 


As  in  a  previous  integration,  we  use  0 oo  as  the  upper  limit  of 
integration  instead  of  the  exact  result.  In  this  case,  we  are 
further  substantiated  in  this  approximation  by  the  fact  that 
near  ©oo~n the  contributions  from  the  above  integrand  to  C&N 
are  near  zero  due  to  the  cos  (p#  factor  in  the  above. 

:/ith  the  previously  calculated  values,  the  shock  shape  can 
now  be  written  as 


Rj«  =  2 


,Go°  ©w  J/V 

(  J0CosY,'%zsl'09+CDZ9*'*,P*  + 

-~i)  Jo  2Sto  / 


^Qrc.  -l5arc+l) 


Also  from  (15.6)  we  have 


Res  -  Rs 


6oo 


Mz. 


r2 -  2£}m _ _ _ 

(feyg-jr^-p  reoo - 

^Ctf*-iXirV  +  ijj  cos^rrc-L 

0  "2.  ©cD 


-4-  1 


V;e  are  now  in  a  Position  to  calculate  layer  properties  based 
on  these  results.  For  the  external  flow  we  have 


po  - 


kith 


Rse  “  C  X 


(#f  *  '-S?  !■  ( 1  -  (ft)-) 

*  p*  -  fe*  ( 1  *  '«*'*) 


This  distrioution  will  be  valid  throughout  the  external  laye: 
i he  density  distribution  is  given  by 
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Po  -  (i  +  ^*/?>)  %/2.v 

The  result  is  not  strictly  valid  since  the  density  is  not 
zero  at  the  contact  surface. 

VJhen  referenced  to  values  at  the  shock,  the  results  become 

P=/P=  =  Va  U  +  HV/P) 

Pe/fe.  ~  Vs.  (HV'+'Xl  +  M4/<)>)  (16.2) 

Te/fe  =  HVy*  **  (16.3) 

For  the  internal  flow,  we  have  the  following  results: 

P«<.  =  C*/ 8x  (16.4) 

foL  -  C’/8X  (16.5) 

Voi  -  \fii/C  £'&/zx,/t  =  M  (16-6) 

When  referenced  to  values  at  the  shock,  the  above  become 


p;/p. 

(16:7) 

?l/p* 

—  — 

'  £g?.iia 

€e 

(16. g) 

®  f^ooc. 

(16.9) 

Vt/0. 

-  % 

C Z  U.OO*. 


**  Not  valid  at  contact  surface 
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The  first  order  solution  obtained  thusfar  from  the  hypersonic 
small  disturbance  form  of  the  IJewton— Busemann  approximation, 
is  adequate  to  describe  layer  thicknesses,  pressures,  and  den¬ 
sities  in  the  far  field  region  near  the  contact  surface.  How¬ 
ever,  due  to  the  assumption  of  &V»i**lthe  temperature  to  the 
first  approximation  is  in  considerable  error  for 
This  is  because  the  expression  for  temperature,^ 
gives  the  unrealistic  result  that  for  all  X/^..  5X^las  *-*1 
This  result  is  acceptable  for  small  ,  ,  where  p/p^  is  of 
order  one,  however,  for  V^1  ^is  small  and  even  for 

near  one  3L  <<1. 

Ts 

To  obtain  a  more  accurate  value  for  temperature  near  the  con¬ 
tact  surface,  the  unexpanded  forms  of  the  variables  in  terms 
of  &  will  be  retained.  Referring  to  Cheng  and  Kirsch3\  the 
unexpanded  for  S”  results  for  the  entropy  wake  are 


R«a.  -  ft*.  =  ^/p'**  (17tl) 

fe  -  (17.2) 

The  retention  of  the  exponent,  'fa ,  in  0.7.1),  (17. 2)  is  justified 

because  the  relative  errors  are  smaller  than  any  integral  po- 

34 

wer  of  .  Applying  the  same  reasoning  to  the  internal  la¬ 

yer  expressions,  (14.3),  (14.4),  we  obtain 

fi  =  (EiL^fi*  (w-3) 
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(17.4) 


Matching  (17.1)  to  the  near  shock  layer  results  (6.1)  and 
then  combining  with  (17.4)  we  obtain 

P^'Cr-s*.  -  R«)  =  <r/a. 

This  is  the  same  form  as  expression  (10.1)  with  po  replaced 

'/jfc  _ 

by  pe  .  To  the  accuracy  needed  for  1  it  is  ade¬ 

quate  to  use  the  formal  result  (10.1)  which  utilises  p0  in- 

,/&£m 

stead  of  pe  .  However,  for  the  distribution  of  temperature, 
we  use 


JL  ^ 
Ts 


(17.5) 


for  both  the  external  near  the  contact  surface  and  internal 


flow.  This  expression  will  give  accurate  values  of 


for 


small  P/P  s  as  well  as  for  P/ps  near  one.  To  be  consistent, 

(17.6) 


we  also  use  the  density  expression 

*  =  (*f 

instead  of  J3  =  £*s  P/p.s 


Making  use  of  (17.5)  and  (17.6)  and  the  expression  for  pres¬ 
sure  (16.1)  the  temperature  and  density  distributions  near 
the  contact  surface  for  the  external  flow  and  within  the  in¬ 
ternal  layer  are 

t  =  T* 


f-dJ) 

\8X  ps/ 

f  = 

\8*p  */ 
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The  above  expressions,  as  well  as  those  previously  found,  are 
used  to  predict  geometry  and  flow  properties  in  the  far  field 
region  for  typical  systems  of  interest.. 
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DISCUSSION  OF  RESULTS 


Examining  the  far  field  analytic  solution,  we  find  that  an 
experimental  simulation  of  an  actual  flow  requires  matching 

*  ©»>  n  and  Rsto  ,  the  same  scaling  parameters  used 
in  the. bow  region  simulation.  If,  in  addition  to  the  above, 

^°J  / Pooe  is  matched  then  the  bow,  corner  and  far  field  geo¬ 
metry  will  be  simulated.  For  a  given  system  with  fixed  am¬ 
bient  and  exhaust  gas  composition  (i.e.,^e,^y  fixed)  and  noz¬ 
zle  conditions  (i.e.,©QD,n  )  the  flow  geometry  will  scale  as 
does  the  bow  region  with  RSt'0  or  TfV*.  if  ri®eis  also  as¬ 
sumed  to  be  fixed,  then  the  bow,  corner  and  far  field  regions 
will  all  scale  with  Rsu»  .  Calculations  of  far  field  flow 
geometry  and  properties  have  been  carried  out  and  the  results 
appear  in  figures  (21)  to  (24).  As  in  the  bow  analysis,  the 
far  field  geometry  is  sensitive  to  the  exponent  used  in  the 
exhaust  plume  model.  From  figures  (21)  and  (22)  it  can  be  ob¬ 
served  that  the  internal  layer  occupies  a  large  fraction  of  the 
shock  layer  flow.  This  is  due  to  the  greater  amount  of  str¬ 
ongly  shocked,  and  hence  lower  density,  gas  in  the  internal 
flow  compared  to  that  in  the  external  flow. 

As  shown  in  figures  (21)  and  (22),  increasing  Boo  or  re¬ 

sults  in  the  moving  outward  of  the  shock  and  contact  surfaces. 
This  behavior  is  consistent  with  the  results  from  the  bow  re¬ 
gion.  Also,  near  the  contact  surface,  temperature,  pressure, 
and  density for  the  internal  and  external  flow  increase  for 
increases  in  Go©  or 


T 


The  effect  of  the  internal  layer  flow  on  the  external  flow 
pressure,  density,  and  temperature  is  to  maintain  them  at  a 
higher  level  than  that  created  by  just  the  bow  interaction. 
This  can  be  seen  from  expression  15.4  in  the  far  field  analy- 
sis  section  where  the  coefficient  for  the  pressure  decay 
expression  consists  of  the  normal  bow  interaction  term  plus 
the  effect  of  inner  layer  thickness. 

In  figures  (23)  and  (24)  the  detailed  property  distributions 
across  the  layers  are  given  for  the  same  "cold"  jet  and  actual 
jet  cases  that  we re  analyzed  in  the  bow  section. 

Consistent  with  the  bow  region  results,  we  find  that  the  act¬ 
ual  jet  case  shown  in  figure  (24)  has  a  hot,  low  density  ex¬ 
ternal  layer,  which  flows  over  a  cooler  and  more  dense  inner 
layer.  If  allowed,  considerable  heat  transfer  and  mixing  bet¬ 
ween  the  external  and  internal  layers  would  occur.  The  "cold" 
jet  case  does  .not  exhibit  as  sharp  a  change  in  flow  proper¬ 
ties  across  the  contact  surface  as  that  of  the  actual  case. 


This  indicates  that  "cold"  jet  experiments  must  be  interpreted 
carefully  in  light  of  the  greater  mixing  and  heat  trans¬ 
fer  effects  for  the  actual  system  versus  the  "cold"  jet  simu¬ 
lation. 
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CONCLUSIONS 


Analytical  techniques  for  predicting  flow  properties  in  the  bow 
and  far  field  regions  of  an  opposed  hypersonic  plume  in  a  hy¬ 
personic  stream  have  been  developed.  It  has  been  shown  that 
they  are  valid  for  a  wide  range  of  altitude  and  .jet  thrusts 
for  typical  systems.  The  influence  of  the  corner  region  flow 
on  the  bow  and  far  field  flow  predictions  has  been  shown  to  be 
negligible.  Consequently,  only  a  qualitative  outline  of  the 
flov;  processes  in  this  region  are  given. 

Some  of  the  major  results  of  the  analyses  and  calculations 
are : 

1.  Experimental  simulation  depends  on  the  matching  of 

i  ^  )  R*w>  )  /  poo e 

2.  For  a  given  system  with  fixed  ambient  and  exhaust 
composition,  external  Ilach  number  and  nozzle  con¬ 
ditions  (i.e.,  fixed  #e,Vj;0eD; n;M®e)  the  entire  flow 
geometry  will  scale  with  Rsi©or  Tf  , 

3.  This  analysis  confirms  the  good  accuracy  of, the 
Newtonian  impact  analysis  when  applied  to  the  axi- 
symmetric  bow  region,  and  suggests  that  the  extension 
of  this  simple  technique  to  predict  asymmetric  bow 
geometry  might  be  successful. 

4.  The  bow  and  far  field  solutions  are  sensitive  to  the 
exhaust  plume  model  exponent  used  and  consequently 
for  accurate  predictions  of  flow  properties  an  accu¬ 
rate  plume  model  must  be  utilized. 
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5.  Penetration  of  the  intercepting  shock  layer  flow  into 
the  internal  layer  is  considerable.  However,  the 
influence  on  the  bow  and  far  field  properties  is 
small. 

6.  Good  agreement  is  found  between  a  calculated  bow  re¬ 
gion  geometry  and  an  experimental  result. 

7.  Calculations  show  that  mixing  and  heat  transfer  ef¬ 
fects  between  the  hot  external  and  cooler  internal 
layers  are  more  pronounced  in  the  actual  case  than  in 
the  experimental  "cold”  jet  simulation. 

8.  Increasing  the  primary  system  design  parameters  0oo 
and  *  with  all  other  parameters  fixed  results  in 

a.  The  thickening  of  internal  and  external  layers 
in  the  bow  and  far  field  regions 

b.  An  increase  in  radius  of  curvature  of  the  con¬ 
tact  surface  and  thereby  an  increase  in  bluntness 
for  the  bow  region 

c.  An  increase  in  the  angle  between  the  far  field 
external  shock  and  free  stream  direction 

d.  A  higher  level  of  pressure  density  and  tem¬ 
perature  both  across  and  along  the  internal  and 
external  layer  flows. 


0(0  is  related  to  nozzle  exit  Mach  number  and  angle  and 
exhaust  gas  composition.  For  a  fixed  composition,  increasing 
nozzle  exit  single  or  decreasing  exit  Idach  number  gives  in¬ 
creases  in  ©oo  .  Vj  is  related  to  exhaust  gas  composition 
which  in  general  will  decrease  as  the  degrees  of  freedom  or 
the  complexity  of  the  exhaust  molecules  increase. 
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9.  The  effect  of  the  internal  layer  on  the  far  field 
region  properties  is  to  maintain  the  pressure ,  den¬ 
sity*  ani  temperature  of  the  flow  at  a  higher  level 
than  that  which  would  be  created  by  the  bow  region 
interaction  with  the  far  field  internal  layer  flow 
absent.  i 

For  high  Reynolds'  number  flows  the  present  technique  may  be 
extended  to  treat  mixing  effects  along  the  contact  surface 
by  the  use  of  boundary  layer  methods.  Also,  equilibrium  chem¬ 
istry  effects  may  be  easily  incorporated  into  the  present  mo¬ 
del  through  use  of  Mollier  charts. 

Further  extension  of  the  bow  layer  technique  to  include  vis¬ 
cous  chemically  reacting  and  merged  layer  effects  requires 
considerably  more  effort  than  the  above  extensions.  In  this 
regard,  an  advantage  of  the  present  technique  over  the  numer¬ 
ical  technique  of  reference  (7)  is  that  a  complete  solution  of 
the  bow  region  requires  approximately  13  seconds  on  an  IBII 
360-65,  whereas  the  numerical  technique  requires  120  seconds 
on  the  much  faster  CDC  6600.  This  economy  may  prove  to  be 
significant  in  extending  the  above  techniques  to  be  able  to 
predict  viscous  chemically  reacting  or  merged  layer  bow  pro¬ 
perties  where  an  order  of  magnitude  increase  in  computational 
times  is  expected.  This  makes  the  technique  of  reference  (7) 
uneconomical  for  use  in  parametric  calculations  whereas  the 
present  technique  would  still  fall  within  the  practical  time 
limitation  for  these  calculations. 
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Appendix  (A) 

Bow  Initializing  Scheme 

To  integrate  the  bow  region  equations  for  points  away  from 
the  axis,  we  must  first  obtain  the  initial  radii  of  curva¬ 
tures  and  positions  of  the  shocks  at  the  axis.  Writing  equa. 
tions  (28.1),  (28.2),  (58. 1)  and  (58.2)  in  the  form 

^  l  3  A  CaO  +  B(q  c)  <5>  t  ( 1 ) 

pL  =  CCa.0  +  b  (fu)  fi*  (2) 

Jjc  »  At<U)  +•  F(Ae')  fe*  (3) 

pc  =  CCaO  +  b'(cie)^ea 

we  develop  expressions  for  and  shock  radii  of  curva¬ 

tures  RsCo/f?s,  .  The  geometry  which  relates  the  internal 
flow  parameters  to  those  of  the  external  flow  is  illustrated 
in  schematic  (5) 


- «■'  - - 
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In  the  above  illustration 


-  Rsco(i  +  A  fee)  +(pfa;))§i)  (5) 

Re  -  Rs*  Cl  +  ACQt')-  BCqc^9e)  (6) 


Referring  to  schematic  (5)  we  can  write  for 
Angle  matching 

fifl  -  +  2.(d«:  4  B'-2AC*A)$t/(l  +  A0  ■+ 

2/<ie  V  aae^A  -  S)fe/tL-S) 

Pressure  Matching 

9e/^c  -(^y 


(7) 

(3) 


V/e  now  have  two  equations  with  three  unknowns  9e/^t*  > 
Consequently,  we  must  relax  one  of  the  parameters  a;,  Qe 
in  order  to  obtain  a  unique  solution.  Choosing  =  0  the 
expressions  for  A  (  0.0,  B  (  Q.C) ,  C  (  <U)  and  D  (  Oj)  simplify 
considerably. 


Substituting  (8)  into  (7)  we  obtain  Qe  which  also  gives  fo/fl. 
From  these  results  the  ratio  of  shock  radii  of  curvatures  at 
the  axis,  R.Si'0/f^St  ,  and  the  relative  position  of  their  cen¬ 
ters  must  be  obtained.  Referring  once  again  to  schematic  (5) 
we  have  from  geometry 

s m  frr-  ^ i )  =i  si>-i  =  s  io  (<? /  -  <Pe3 

Re  R.  ’ 

For  small  values  of  the  angles  we  can  write 

-  Rsg-  C  1  ~  A  ~h  (fle  ~  EOpe»*~) 

Pe  ^(TTMWr 
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For  the  initial  point,  can  be  obtained. 

The  distance  between  their  centers  is  then  given  by 


Having  obtained  the  ratio  of  the  radii  of  curvatures  and 
their  relative  positions  at  the  axis,  we  then  use  these  quan 
tities  as  input,  data  to  the  computer  program,  which  predict 
bow  properties  for  regions  away  from  the  axis  of  symmetry. 


SED  ON  NORMAL  SHOCK  RELAXATION 
R  EXTERNAL  FLOW  (14) 


FIG.  2  FLOW  REGIMES 
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FIG.  3  UNDISTURBED  PLUME  MODEL  EXPONENT  AS  A 
FUNCTION  OF  ANGLE  AND  DISTANCE 


FIG.  4  EXPERIMENTALLY  DETERMINED  SHOCK  ANO 
CONTACT  SURFACE  LOCATIONS  REFERREO 
TO  R, 
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PREDICTED  CONTACT  SURFACE  SHAPES 


COMPARISON  OF  SHOCK  AND  CONTACT  SURFACE  POSITIONS 
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FIG.  7  COMPARISON  OF  CONTACT  SURFACE  PRESSURE  DISTRIBUTIONS 
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VARIATION  OF  BOW  GEOMETRY  WITH 
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FIG.  12  VARIATION  OF  CONTACT  SURFACE  PRESSURE  WITH 


VARIATION  OF  BOW  GEOMETRY  WITH 


FIG.  14  VARIATION  OF  CONTACT  SURFACE  PRESSURE  WITH 
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ISOBARS  AND  STREAMLINES  FOR  BOW  REGION 


16  DENSITY  AND  VELOCITY  DISTRIBUTIONS 


LIM ITING 


FIG.  19  CORNER  REGION 


MACH  NUMBER 


FIELD  SHOCK  AND  CONTACT  SURFACE  POSITION 


CONTACT  SURFACE  POSITION 


FIG. 23  FAR  FIELD  PROPERTIES  VS.  STREAM 
FUNCTION 


EXTERNAL  LAYER 


